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ABSTRACT 

Using cosmological hydrodynamic simulations that dynamically incorporate enriched 
galactic outflows together with analytical modeling, we study the origin of the stellar 
mass-gas-phase metallicity relation (MZR). We find that metallicities are driven by 
an equilibrium between the rate of enrichment owing to star formation and the rate of 
dilution owing to infall of unenriched gas. This equilibrium is in turn governed by the 
outflow strength. As such, the MZR provides valuable insights and strong constraints 
on galactic outflow properties across cosmic time. We compare t hree outflow models : 
No outflows, a "constant wind" model that emulates the popular fPekel fc Silkl (|1986D 
scenario, and a "momentum-driv en wind" model that best r eproduces z > 2 inter- 
galactic medium metallicity data (jOppenheimer fc Davel[2006h . Only the momentum- 
driven wind scaling simulation is able to reproduce the observed z ^ 2 MZR's slope, 
amplitude, and scatter. In order to understand why, we construct a one-zone chem- 
ical evolution model guided by simulations. This model shows that the MZR in our 
outflow simulations can be understood in terms of three parameters: (1) The equilib- 
rium metallicity Zg^eq = ?/-^SFR/AfACC (where j/=net yield), reflecting the enrichment 
balance between star formation rate Msfr and gas accretion rate Macc; (2) the di- 
lution time t(i — Mg/Mf^cQ, representing the timescale for a galaxy to return to 
Zg.cq after a metallicity-perturbing interaction; and (3) the blowout mass Mbiowout, 
which is the galaxy stellar mass above which winds can escape its halo. Without out- 
flows, galaxy metallicities exceed observations by ~ x2 — 3, although the slope of the 
MZR is roughly correct owing to greater star formation efRciencies in larger galaxies. 
When outflows with mass loading factor jyw are present, galaxies below Mbiowout obey 
Zg,eq ~ + w)i while above Mbiowout, Zg^cq ^ U- Our constant wind model has 
Mbiowout ^ 10^*^ M0, which yields a sharp upturn in the MZR above this scale and a 
flat MZR with large scatter below it, in strong disagreement with observations. Our 
momentum-driven wind model naturally reproduces the observed Zg oc M^-^ because 

Zg,eq oc r]^ ex. mI^^ when t/w ^ 1 (i.e. at low masse s) . The flattening of the MZR 
at M* > 10^°-^Mq observed by iTremonti et all (|2004l ) is reflective of the mass scale 
where rp^ ~ 1, rather than a characteristic outflow speed; in fact, the outflow speed 
plays little role in the MZR except through Mbiowout- The tight observed MZR scatter 
is ensured when td < dynamical time, which is only satisifled at all masses in our 
momentum-driven wind model. We also discuss secondary effects on the MZR, such 
as baryonic stripping from neighboring galaxies' outflows. 
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1 INTRODUCTION 

In the reigning hierarchical model of galaxy formation, cool- 
ing times for moderate overdensity gas at high redshifts 
{z > 3) are shorter than a Hubble time, and protogalax- 
ies accrete gas from the intergalac tic medium (IGM) a t 
rates that increase with time (e.g., [Birnboim et al.|[2007l ). 



Gas forms into stars a t rates that main ly track the ris- 
ing gas accretion rates (|Keres et al.l 120051 ): these rates also 
scale with g alaxy mass and are o c casionally boosted b y 
mergers (e.g. ^ Somerville et al. I I2OOII : IStraughn et"alll2006l ). 
Each generation of stars in turn enriches the surrounding 
gas with heavy elements. The young galaxies also drive 
much of their enriched gas out via galactic superwinds, lac- 
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ing the IGM with he avy elements (jAdelberger et al]l2005bl . 
I2OO3I : ISchaereill2003l) while suppress i ng their own star for- 
matio n ( Springel fc Her nguist 2 003bl : lOppenheimer 
I2OO6I ). As the universe expands further, cooling times rise, 
until z ~ 1 — 2 when the coohng time at moderate overden- 
sities exceeds the Hubble time, and gas accretion rates and 
star formation rates begin to decline. 

At any given epoch the accumulated history of star for- 
mation, inflows, and outflows affects a galaxy's mass and 
its metallicity. Hence one expects these quantities (or their 
proxies) to be correlated in some way, and furthermore for 
this correlation to encode information about the physical 
processes that govern galaxy formation. In this paper, we in- 
vestigate what constraints may be placed on such processes, 
particularly outflow processes that are currently the most 
poorly understood, based on the observed mass-metallicity 
rela tion (MZR) of galaxies. 

iMcClure fc van den Berghl ()l968l l were the first to ob- 
serve a correlation bet ween the lumino s ity an d metallicity of 
elliptical galaxies, and iLeaueux et al.l l| 19791 ) were the first 
to show that total mass (or, equivalently, rotational veloc- 
ity) correlates with metallicity for irregular and blue com- 
pact galajcies. The question of which relationship is more 
fundamental dZaritskv et al."l994': 'G arnett et al.lll997l ) was 
resolved when Trcmonti ot al. (2004) showed, for a sample 
of « 53, 000 star-forming galaxies from the Sloan Digital 
Sky Survey, that the MZR possesses much less intrinsic scat- 
ter thBB_th«jAimino relation. The trends seen 
bv lTremonti eUD |200i) have recently been shown to ex - 
tend unbroken to much lower masses by iLee et all ||2006D . 
confirming the idea (jGarnettl I2OO2I ) that a single mecha- 
nism may govern galaxies' metallicities across five decades in 
stellar mass. Finally, the MZR is observed to evolve slowly 
such that galaxies of a given stellar mass arc only moder- 
ately r norc enriched today as compared to in the e arly Uni- 
verse llShaplev et al.l l2005l : ISavaglio et all l2005l : lErb et al.l 
l2006l : lBerger et al.ll2006h 

A number of mechanisms have been proposed to ex- 
plain the observed trends. Building upon a n ide a origi- 
nally int roduced by j Mathews fc Bakei] lll97ll) and iLarsonl 
l|l974h . iDekel fc Silkl l|l986l ) and iDekel fc Wool (|2003l ) 
showed that supernova feedback energy could give rise to 
a range of observed trends in low mass galaxies including 
the MZR, if the supernova energy injected into galaxies' in- 
terstellar media is proportional to its stellar mass (as one 
might naively expect). In this model, there is a character- 
istic halo circular velocity of ~ 100 km s~^ above which 
galaxies retain their gas, and below which gas removal be- 
comes progressively more efficient. 

More recent investigations have attempted to pu t the 
MZR in a hierarchical context. iDe Lucia et al.l (|200J) used 
semianalytical models to suggest that such winds can be 
tuned to reproduce the 2 « MZR irres pective of what 
the o utflo ws do upon l eaving the galaxy. iDe Rossi et al.l 
l|2006t ) and lTassis et al.l l|2006l ) used cosmological hydrody- 
namic simulations without strong supernova feedback to 
obt ain ro ugh agreement with the iTremonti et al.l (|2004l ) 
and iDeke l & Woo (200i5i) MZRs, respectively. Both works 
cite the possible role of a varying star formation efficiency 
with stellar mass, while iTassis et al.l (|2006l l also suggest 
that the observed low effective yields in low-mass galaxies 
could result from mixing processes that transport metals 



to galaxies' outer d isks, where they are difficult to observe. 
iBrooks et aP l|2006l ) used cosmological hydrodynamic sim- 
ulations with a treatment for pressure-driven outfiows to 
argue that mass loss does not directly suppress the metal- 
licities of low mass galaxies, and instead argued that su- 
pernova feedback leads to low star formation efficiencies in 
low-mass galaxies, which in turn lead to low metallicities; it 
is particularly encouraging that their mod el reproduces the 
observ ed trends both at 2: ~ 2 and at 2 ^ O. lKobavashi et ahl 
iiooj) used a similar model that also incorporated a treat- 
ment for hypernova feedback and obtained qualitative agree- 
ment with the same observations. However, they concluded 
that their MZR is primarily driven by a tendency for low- 
mass galaxies to eject relatively more material in outfiows, 
an i dea that our results support. A more speculative idea 
from lKoppen et al] (I2OO6I ) suggests that the observed trends 
at low redshift can be reproduced by postulating that the 
stellar initial mass function (IMF) is more top-heavy in 
galaxies with higher star formation rates, thereby producing 
higher metal yields. 

In short, the galaxy MZR has been speculated as arising 
from variation in mass loss, star formation efficiency, and/or 
yield with galaxy stellar mass. A key aim of this present work 
is to distinguish between these alternatives, and determine 
the key driver that sets the MZR. 

A feature that has garnered much attention is the ap- 
parent flattening of the MZR for M* > 10^°-^ Mq. This char- 
acteristic mass also seems to divide galaxy properties in gen- 
eral, such as blue from red and high surface brightness from 
low (Kauffmann ot al. 2004). Dckcl fc Woo (200 ^) note that 
their predicted characteristic halo velocity ( Dekel fc Silkl 
198611 is in broad agree ment with this mass. B uilding on this. 



Tremonti et~ai1 l|2004l ) and jGarnettl l|2002h proposed that 



winds may not be effective at driving me tals out of galax- 
ies above this mass range (see however iDalcantoru i2006l ). 
One physical model that could give rise to this behavior is 
the "constant wind" scenario. In this model, the low effec- 
tive yields observed in small galaxies result from winds that 
have roughly constant velocities at all masses, and hence are 
progressively more effective at driving out metals from the 
shallower potential wells of smaller galaxies. Consequently, 
the escape velocity at the observed MZR turnover should 
be an indicator of the characteristic wind speed. At present, 
this scenario of mass loss variation is probably the most 
widely accepted explanation for the MZR's shape. 

In this work we show that such a constant wind sce- 
nario, when incorporated into a fully three-dimensional hi- 
erarchical structure formation simulation, produces an MZR 
shape that is in poor agreement with observations, for rea- 
sons that can be understood from straightforward physical 
arguments. This scenario has been also called into question 
recently by obs ervational analyses of dwarf galaxy metal- 
licities. iLee et al. (2006) note that such energetic speeds 
from small galaxies should produce a much greater scatter 
in metallicities than observed in their sample of dwarf irreg- 
ulars. They instead propose that "a less energetic form of 
metal-enhanced mass loss than blowouts could explain the 
small scatter." ,Dalcantonl l|2006l) used a rigorous treatment 
of the effects of outflows to show that outfiows alone cannot 
account for the low observed effective yields of dwarf galax- 
ies, unless the winds are substantially enriched relative to 
ISM gas. She proposed instead that the low gas surface den- 
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sities of galaxies with circular velocities below 120 km s~^ 
lead to low star formation rates (SFRs) so that their effec- 
tive yields "recover" from outflows relatively slowly. Hence 
she suggests star formation efliciency variations are the key 
driver of the MZR. The results we present here are generally 
in agreement with this conclusion. 

In this paper we use cosmological hydrodynamic simu- 
lations and simple analytical models to investigate the ori- 
gin of the mass-metallicity relation. Our simulations em- 
ploy parameterized outflows from star-forming galaxies that 
drive metals into the IGM, and hence we directly track the 
growth of galaxy metallicity along with its mass within a 
full three -dimensional hierarchic a l stru cture formation sce- 
nario. In lOppenheimer fc Pavel l|2006l ) we introduced our 
outflow models, and showed in particular that models in 
which the outflow velocity scaled linearly with galaxy cir- 
cular velocity while the mass loading factor (i.e. the rate 
of mass ejection relative to the star formation rate) scaled 
inversely with it were remarkably successful at enriching 
the IGM to observed levels at z ~ 2 — 6. These scalings 
arise naturally for radiation o r momentum-driven winds (e.g. 
iMurrav. Quatert. fc Thompso n 2005), though for our pur- 
poses the important aspect is the scaling relations them- 
selves and not the physical mechanism responsible. Interest- 
ingly, recent observations of local galactic outflows indicate 
momentum-driven wind scalings (|Martinll200"5l : lRupke et al.l 
l2005h . providing an intriguing connection between rare lo- 
cal outflows and the more ubiquitous and generally stronger 
outflows at z > 2. 

In our preliminary study of the MZR (jPave et al.l 
l2006bl ). we co mpared var i ous o utflow models with the z ^ 2 
MZR seen bv lErb et al.l (|2006l ) and examined its evolution 
from 2 = 6^2. W e found that the "constant win d" model 
as implemented by ISpringel fc Hernguistl l|2003bD leads to 
poor agreement with observations aX z ~ 2, while our 
momentum-driven wind scalings naturally reproduce the 
slope and amplitude of the observed relation. While this 
work provided independent support for the outflow model 
concurrently favored by IGM metallicity data, it did little to 
address the fundamental question of what physical processes 
govern the MZR's slope, amplitude, scatter, and evolution. 

In order to improve our understanding of the relation- 
ship between outflows and the MZR, we follow a trajectory 
that encompasses three basic goals: 

(1) Show that our numerical model reproduces observations 
at z 2. In this step we describe our simulations (Section[21) 
and discuss the mass scales that our outflow models intro- 
duce. Next we compare the simulated MZRs assuming differ- 
ent outflow models with the observed MZR at z = 2 (Sec- 
tion |3]). Our finding that momentum-driven outflows pro- 
duce the best agreement with observations motivates a more 
detailed investigation into the origin of the MZR within our 
simulations. 

(2) "Boil down" our numerical model to a set of key pro- 
cesses and combine them in an analytically tractable way. 
We begin this step by introducing a simple model that cap- 
tures the main processes that impact the growth of galaxies' 
stellar masses and metallicities (Section |4]). Next, we inves- 
tigate how gas accretion and star formation rates vary with 
mass and time in each outflow scenario in order to treat 
them accurately in our analytical model. We continue this 
discussion in Section [5] with a detailed investigation into the 
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''Box length of cubic volume, in comoving 
ft-^Mpc. 

''Equivalent Plummer gravitational softening 
length, in comoving /i~^kpc. 

'^AU masses quoted in units of 10^ Mq. 

'^Minimum resolved galaxy stellar mass. 

°First number for nw, cw runs; second for vzw. 

time-integrated effects of outflows on our simulated galax- 
ies at 2: = 2. In Section [6] we trace the observable trends 
at 2 = 2 back in time in order to understand how galaxies 
evolve through the MZR. In Section [6.31 we verify that our 
analytical model reproduces this evolution reasonably well, 
which suggests that our analytical model accounts for the 
important processes that drive the observable MZR. 
(3) Use the analytical model to determine what drives the 
observable MZR in the presence of outflows. In Section [7] we 
use our analytical model to determine what drives the form 
of the observable MZR; the reader may wish to skip directly 
to this Section for a relatively self-contained explanation of 
the origin of the MZR. Here we show why our momentum- 
driven wind model is successful at reproducing the z « 2 
MZR, and why other models are not. More generally, we 
show how the observed MZR's slope, amplitude, and scatter 
down to low masses provide stringent constraints on outflow 
models. 

Finally, in Section [8] we present our conclusions. 



2 SIMULATIONS 

2.1 Simulations and Sample Definition 

We employ th e parallel cosmological ga laxy formation 
code Gadget-2 (|Springel fc Hernguistl 120021 ') in this study. 
This code uses an entropy-conservative formulation of 
smoothed particle hydrodynamics (SPH) along with a 
tree-particle-mesh algorithm for handling gravity. Heating 
is inclu ded via a spatially un iform photoionizing back- 
ground (jHaardt fc Madaull200lh . Gas particles undergo ra- 
diative cooling under the assumption of ionization equilib- 
rium, where we account for metal-line cooling usi ng the colli- 
si onal ionization equilibrium tables of Sutherlan d fc Dopital 
(| 19931 ). The metal cooling function is interpolated to the gas 
metallic ity as tracked self-consistent ly by Gadget-2 (for de- 
tails see lOppenheimer fc Dave|[2006h . Stars are formed from 
dense gas via a subresolution multi-phase model that tracks 
conde nsation and evaporation in the interstellar medium fol- 
lowing [McKee^^fc^Ostrike^ 1 19771 ). The model is tuned via a 
sing le paramete r , the s tar formation t imescale, to reproduce 
the iKennicuttI l|l998l ) relation; see ISpringel fc Hernguistl 
(|2003al ) for details. Star-forming gas continually self-enriches 
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under an instantaneous recycling approximation. The ef- 
fects of Type la supernovae are neglected; this should not 
affect comparisons with observed Oxygen abundances. In 
reality Type II supernovae occur with a time delay of 10- 
30 Myr, longer than our typical timestep of ~a few Myr, 
hence the assumption of instantaneous feedback is inappro- 
priate if gas accretion or star formation is believed to vary 
on shorter timescales. Howev er, our simulations (and ob- 
servations; iNoeske et aDl2007l l indicate that star formation 
occurs in a predominantly smooth fashion, hence instanta- 
neous feedback is unlikely to introduce significant errors. De- 
layed feedback from low-mass stars is also neglected. When 
star particles are spawned (in two stages, each with half the 
original gas particle's mass), they inherit the metallicity of 
the parent gas particle and from then on cannot be further 
enriched. 

We model galactic outfiows using a Monte Carlo ap- 
proach. First we define the wind model by two parameters: 
a wind speed Vw and a mass loading factor r/w, which is 
the ratio of the outfiow mass rate to the star formation rate. 
These parameters can be chosen to be constant or scale with 
galaxy properties. During the simulation run, for each star- 
forming particle we compute a probability that it enters into 
an outflow based on its star formation rate and r/w, and use 
a random number to decide whether that particle will enter 
into an outflow. If so, we kick it with a velocity of Vw in the 
direction of vxa, which would be purely unipolar for a thin 
disk but more typically has a large opening angle of ~ 45°. 
The hydrodynamic forces are turned off for that wind parti- 
cle until it reaches a density of one-tenth the critical density 
for star formation, or else it travels for a time greater than 
(20 kpc/h)/Vw. The outfiow particle carries its own metals 
out of the galaxy; it is not preferentially enriched. 

In order to explore the effects of superwind feedback 
on the MZR, we concentrate on three outflow schemes: A 
no-wind (nw) model, a "constant wind" (cw) model where 
Vw = 484 kms~^ and riw = 2 (this is the scheme 
used in the runs of ISpringel fc Hernguisd l2003bl ). and a 
"momentum-driven wind" (vzw) model where Vw cx; a 
and rj = (300 'kms~^)/a, where the velocity disper- 
sion a is estimated from the local gravitational potential. 
The exact scali ngs are taken from the momentum-d riven 
wind model of Murray. Quatert. fc ThompsonI (120051 ): see 
lOppenheimer fc Davel li20m for details. The vzw model also 
gives a velocity boost in low-metallicity systems, based on 
the arguments that more UV photons are produced per unit 
stellar mass at lower metallicities (specifically, we employ 
eqn. 1 of lSchaereil2"oo|V, and that it is these UV photons tha t 
are driving the wind (^ Murrav. Quatert. fc Thompsonll2005h . 
Using larger simulations evolved to low redshift, we have 
verified that both models bro adly reproduce the star for- 
mation history of the universe ijOppenheimer fc Davellioo^ . 
Figure 4); this constrains the choice of parameter values. 

We run a total of six simulations: two different vol- 
umes, each with our three different superwind schemes, as 
detailed in Table [1] All runs a ssume a WMAP-concordant 
cosmology fepergel et al]|2003l ) having Q = 0.3, A = 0.7, 
Ho = 70 kms"^ Mpc"\ erg = 0.9, and fit = 0.04. Each run 
has 256^ dark matter and 256'^ gas particles, evolved from 
well in the linear regime. 

We identify galaxies using Spline Kernel Interpolative 
DENMAX and dark matter halos using the spherical over- 



density algorithm (see iKeres et al.l |2005| . for full descrip- 
tions). In previous papers we have shown that imposing 
a mass resolution cut of 64 star particles leads to a con- 
verged sample in terms of stellar mass and star for mation 
history (jFinlator et al.ll2006l . [2007l : lDave et"aLll2006al '). In the 
present work we make an even more conservative cut at 128 
star particles in order to study the scatter in our simulated 
trends. Hence our minimum resolved galaxy stellar mass 
ranges from 2.5 x 10* M© in our higher-resolution 16/i~^Mpc 
runs to 2.0 X lO^Afo in our la r ger vol ume S2h ^Mpc runs. 
Just as we did in |d ave et al.l (|2006al ). we will often show 
both runs on the same plot, and the smoothness of trends 
in overlapping mass ranges is to be noted as an indicator of 
numerical resolution convergence. We have also run higher- 
resolution simulations in 8/i~^Mpc volumes and confirmed 
that the trends that we identify continue to lower masses, 
though the galaxies in these runs are unobservably small at 
the redshifts where we will perform comparisons to data. 

2.2 Scales in the Wind Models 

Our wind models introduce several mass scales that are im- 
portant in understanding the behavior of the MZR. The 
reheating scale is the scale below which galaxies produce 
enough supernova energy to unbind all of their gas. The 
blowout scale is the scale below which the wind speed exceeds 
the escape velocity of the halo. Here we compute values for 
these scales in our wind models. 

We first consider the reheating scale in the cw model. 
The virial energy of the baryons in a halo scales with the halo 
mass as i?vir oc M^^^ whereas the feedback energy scales as 
Evind OC M (assuming that the fraction of baryons converted 
to stars /, varies slowly with M). The ratio of these energies 
thus scales as -Bwind/^?vir cx: M~^^^ with the implication that 
the relative importance of wind heating declines as mass 
increases in the cw model. Below the "reheating scale", a 
galaxy's winds produce enough energy to expel all of the 
baryons from the halo; this is analogous to the critical scale 
for supernova-driven mass loss proposed by iDekel & SiM 
(198 6). In the spherical collapse scenario (e.g., iDekel fc Wool 
[2003) it can be shown that this scale corresponds to a stellar 
mass of 
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Here, A(2) ~ 200 is the overdensit y of collapsed struct ures, 
which varies weakly with redshift (jPekel fc Woci|[2003h : ?7w 
is the mass loading factor of the winds; and Vw is the wind 
speed. Note that the circular velocity of the reheating scale 
in the cw model is sim ilar to the critical velocity identified 
bv lDekel fc Sii5 ()l986l l. This is expected because the energy 
injected by winds assuming these p arameters is compara- 
ble to the supernova feedback energy (jSpringel fc Hernguistl 
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l2003bl ). Hence a comparison between our cw model and ob- 
servat ions constitutes a quantitative test of the lDekel fc Silkl 
l|l986l ) scenario. 

In galaxies above the blowout scale, our constant wind 
model outflows should escape from the galaxies' halos pro- 
vided that they couple inefficiently with the ambient halo 
gas. Under reasonable assumptions regarding the depth of 
the galaxy's potential well, the ratio of the escape velocity 
to the halo c ircular veloci ty lies within the range Vasc/V ~ 
f.5-3.5 (e.g.. iMartinlligggl ). In the spherical collapse model, 
the blowout scale in the cw model is therefore given by 
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Both the reheating scale and the blowout scale fall 
within the range that is resolved by our simulations, hence 
we expect to see features in the MZR indicating that the 
two effects grow significantly more effective below this scale 
and decreasingly effective above it. In Section l5.f I we will 
show that this is indeed what happens, although we will ar- 
gue that blowout is much more important than reheating in 
our models. In Section [7] we will use simple physical argu- 
ments to show how the existence of a blowout scale leads to 
predictions that conflict with observations. 

In the vzw model, 77W = co/cr and Vw = fccr, where 
ao is a constant, a is the halo velocity dispersion, and k 
relates the wind velocity and the halo velocity dispersion to 
the ratio of the galax:y luminosity to the galactic E ddington 
luminosity (iMurrav. Quatert. fc Thompson! l2005h : k — 6.7 
on average. Hence by construction, all galaxies are above 
the blowout scale in this model. With these assumptions, 
the injected energy scales as -Ewind oc Alriwcr^ oc a'' cx M*^'^ , 
and the ratio of the heating energy to the virial energy scales 
as E^ind/Evii oc M~^^^ — more shallowly than in the cw 
model. In this case, the reheating scale is given by 
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which is so large that in practice all of our simulated galax- 
ies would expel all their baryons if wind energy coupled ef- 
ficiently with the remaining gas in the galaxies' halos. In 
Figure [7] we will show that this does not happen, and halos 
actually retain much of their baryonic mass far below the re- 
heating scale. Hence in our models, the outflow energy does 
not couple efflciently to the ambient ISM or halo gas, but 
rather tends to blow holes and escape into the IGM. This 
idea is qualitatively consistent with hig h-resolution individ- 
ual g alaxy simulations of outflows (e.g. iMac Low fc Ferraral 
I1999I ). but inconsistent with the assumption of efficient en- 
ergy coupling with ambient gas that is sometimes made in 
analytical blowout models. 
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Fi gure 1. The M ZRs at 2 = 2 co mpared to observations 
bv TErb et aP l|2006l ') (reproduced from IPave et al.l r2006b'l. The 
no-wind case overproduces metals, the constant wind case shows 
too steep a slope above the blowout scale and a large scatter be- 
low it, and the momentum-driven wind scenario fits observations 
quite well. The two clumps of points in each figure correspond to 
the 16 and 32?t~'^Mpc simulation volumes, and are bounded at 
the low-mass end by their galaxy mass resolution limits. 



3 THE Ah - Zg relation AT Z = 2 
3.1 Gas-phase Metallicities 

We begin by demonstrating that the choice of wind model 
heavily impacts the simulation's predictions for the MZR. 
Before doing so, several remarks are in order regarding 
the simulated and obse rved measur e ments. First, through- 
out this work we follow iDave et all (|2006bl ) and define the 
metallicity of each simulated galaxy as the SFR-weighted 
metallicity of its gas particles; this presumably provides 
a fair analogue to metallicities derived from metal emis- 
sion lines as these trace the most actively star-forming 
regions in each galaxy. Second, in order to co mpare our 
simulated metallicites with the measurements of lErb et al.l 
(|2006l ) in Figures [1] and llll we normalize all metallici- 
ties to the net yield. We have calculat ed the net Oxygen 



yield using published Type II SN yields llWoosley fc Weaver] 
I1995I: IChieffi fc Limondl I2OO4I: IPortinari et al.1 Il998h with 
the lSalpeterl (|l955l ) and ichabrieil l|2003D IMFs and find that 



it lies between 0.008 and 0.021 depending on the choice of 
models, IMF, and metallicity. Our simulations assume a to- 



6 Finlator & Dave 



tal metal yield of 0.02 with solar abundance ratios, corre- 
sponding to a net Oxygen yi e ld of 0.0087. We use this value 
to normalize the lErb et aU (l2006l ) measurements for con- 
sistency with our simulations while noting that the choice 
of net yield introduces a factor of ~2 uncertainty in the 
relative normalizations in Figure [T] that is in addition to 
the factor of ~ 2 unce rtainty in the observed metallici- 
ties (IShapley et al ."2005') . Finally, because the stellar masses 
reported in Ikrb et a l. (200^) correspond to the total stellar 
mass formed whereas our simulations report the mass of 
stars that does not im mediately explod e in Type II super- 
novae, we multiply the lErb et al.l (|2006l ) masses by the ratio 
of mass in long- lived stars to t he total mass formed (which 
is 0.802 for the IChabrie"ill2003l IMF) while noting that the 
observed stellar masses are also uncertain at the ~ 2x level 
owing to uncertainty in the IMF. 

Figure [T] shows the MZR that arises i n each wind mode l 
at 2 = 2, compared with observations by lErb et all l|2006t ). 
The fact that outflows affect the MZR is abundantly clear 
from these figures. In the no wind case, galaxies are too en- 
riched at a given stellar mass (this is also seen in the wind- 
free models of iDe Rossi et al.ii2006i ). This indicates that out- 
flows are necessary to expel metals from galaxies, adding to 
the growing body of evidence that outflows have a signifi- 
cant and ubiquitous impact on high-redshift galaxies. Inter- 
estingly, the slope of the MZR in the no-wind case is in fair 
agreement with observations (though slightly too shallow 
in detail), despite the fact that no galaxies are driving out 
any metals. Metal loss by tidal stripping is generally negli- 
gible, as we show in i|5.3l Hence the slope of the MZR does 
not necessarily imply that low-mass galaxies preferentially 
expel metals relative to high-mass ones. Even without out- 
flows, the observed MZR slope (though not its amplitude) 
is broadly reproduced. 

The constant wind (cw) model shows some remarkable 
and interesting trends. Above the blowout scale of IO'^'^Mq, 
the MZR slope is quite steep, steeper than the no- wind case, 
and the scatter is fairly small. Hence relative to the no- 
wind case, the cw model appears to be preferentially ejecting 
metals from low-mass systems, as expected. Unfortunately, 
this takes what was a good agreement with the observed 
slope in the nw case and produces poor agreement in the 
cw case. Below the blowout scale, the scatter increases dra- 
matically; this behavior is seen at all redshifts down from 
2 r; 6 — > and is not a consequence of numerical reso- 
lution, as evidenced by the fact that our higher-resolution 
run joins smoothly onto the lower-resolution one. The large 
scatter in the c w model can be q ualitatively compared with 
the findings of iGeha et al.l (2006), who report that the ob- 
served baryonic Tully-Fisher relation does not show the ex- 
cess scatter below the blowout scale that would be pre dicted 
in blowout scenarios. Additionally, iLee et all (|2006l ) have 
shown that the observed scatter is « 0.1 dex at all mass 
scales at low redshift and argued that this is inconsistent 
with blowout scenarios. We discuss the likely source of this 
excess scatter in Section FT.SI the important point is that this 
scatter is not observed, hence models that introduce a char- 
acteristic blowout scale at observab le masses are unlikely to 
be consistent with iLee et al.l (|2006l ). 

The vzw model's MZR shows good agreement with ob- 
servations in both slope and amplitude. There is a minor off- 
set in amplitude, but there are enough systematic uncertain- 
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Figure 2. Zf versus M, at 2 = 2. Red dots denote galaxies in 
our simulations while blue curves denote the mean trend from the 
gas-phase MZR (Figure[T]l. The data point with error bars denotes 
the inferred stellar mass and metallicity of MS 1512-cB58. Within 
each wind model, the UV-weighted stellar metallicities closely 
track the metallicity of the star-forming gas, as expected. 



ties in met allicity indicators to render this difference insi gnif- 
icant (e.g. lTremonti et al .''2004': 'Ellison fc KewlevlbOOSl ). To 
explain this agreement, _Dav6 ct al. (2006b) proposed that 
galaxies must lose a roughly fixed fraction of their metals 
independent of their mass, because this would preserve the 
slope seen in the nw case while lowering the amplitude. In- 
deed, by comparing the mass of metals produced by each 
galaxy in the vzw run with the mass of metals retained in 
stars and gas, we find the fraction of metals that galaxies 
retain scales very slowly as M^ °^ at 2 = 2 (tj5.3l). It is in- 
triguing to compare this result to iGeha et al.l ( 20061 ). who 
infer that the fraction of baryons lost cannot vary with bary- 
onic mass based on the fact that the slope of the observed 
~ baryonic Tully-Fisher relation lies very close to the 
value that is expected from cosmological simulations that 
do not treat baryons. In the momentum-driven wind sce- 
nario this arises naturally because outflow speeds scale with 
galaxy escape velocities. 



3.2 Stellar Metallicities 

Although the goal of this paper is to understand the re- 
lationship between stellar mass and gas-phase metallicity. 
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recent observational efforts to constrain the evolution of 
the stellar mass-stellar metallicity relation (hereafter, the 
"stellar MZR") motivate a comparison between our sim- 
ulated gas-phase and stellar metallicities. Before we pro- 
ceed, however, a few comments on systematic effects are 
in order. The first observational constraints on the high- 
redshift stellar MZR will focus on rest-frame UV a bsorption 
features (e.g., the 1978 index of iRix et al.l l2004h because 
at 2; J5 2 these features are redshifted into conveniently- 
accessible optical wavelengths. Unfortunately, they are also 
dominated by sta rs with zero-ag e main-sequence masses 
greater than SM© l|Rix et al.ll2004r ) , hence they are only sen- 
sitive to a galaxies' youngest stars. In order to estimate the 
extent to which the observed stellar MZR depends on the 
choice of metallicity indicator, we have measured the stellar 
MZR for each wind model twice: once using mass-weighted 
stellar metallicities averaged over all stars (Z*,aii), and once 
using only stars that are younger than 100 Myr (^,,uv). The 
latter figure approximates a UV-luminosity-weighted stellar 
metallicity. We find that, for all three wind models, ^,,uv 
lies systematically 40-60% higher than .Z,,aii, indicating that 
UV indices systematically overestimate stellar populations' 
mean metallicities even in the absence of any systematics in 
the indices themselves. We also find that, in all three wind 
models, the scatter in Z*,uv at given A/* is 0.1-0.2 dex, 
whereas the scatter in Z*,an is generally ~ 0.05 dex, tighter 
than the scatter in the gas-phase MZR. The increased scat- 
ter in Zt,\jY likely owes partially to our simulations' limited 
mass resolution. On the other hand, it is also likely that the 
short time baseline sampled by UV indices renders them 
more sensitive to the distribution of metallicities among in- 
dividual HII regions within individual galaxies, hence we do 
expect the scatter in the measured MZR to be larger for UV 
indices than for optical indices. In order to facilitate com- 
parison with upcoming measurements at high redshift, we 
only consider Zt,uv throughout the rest of this work. 

The red points in Figure [2] compare the predicted stel- 
lar MZRs a,t z — 2 for our three wind models (red points) 
to the mean gas-phase MZRs (blue curves) as well as to 
current constraints on the z — 2.7276 lens ed galaxy MS 
1512-CB58 (hereafter, "cB58 "lYee et al.lll996l ). For cB58, we 



use the stellar metallicity of iRix et al.l ( 20041) normaliz ed to 
solar yield and the stellar mass of iBaker et al] l|2004h . We 
assume a factor of 2 uncertainty for each estimate. Com- 
paring the stellar and gas-phase MZRs in Figure [5] indicates 
that they are predicted to be quite similar, regardless of 
the outflow scenario. This is expected since young stars' 
metallicities should track the metallicities of their parent 
gas clouds. On the other hand, it is not consistent with ob- 
servations of cB58, for which the most likely stellar metal- 
licity Z,—0.7Zq significantly exceeds the inf erred gas-phase 
metallicity Zg = 0A±0.1Zq |RixeiaL|[200i). However, con- 
sidering that the metallicity offset in cB58 is well within the 
range of systematic uncertainties and that its ISM metallic- 
ity shows excellent agreement with the observed mean gas- 
phase metallicity for its stellar mass at 2; ~ 2 (jErb et al.l 
I2OO6I ). we believe that the offset likely results from system- 
atic offsets in the observational abundance indicators. It will 
be interesting to see whether larger samples of galaxies at 
high redshift show a similar offset. 

We have additionally compared our stellar MZRs to the 
stellar MZR predicted at z = 2 by the hypernova feedback 



model of IKobavashi et al.l (12007 , Figure 20). In the range 
10^~"-^M^. the IKobavashi et al l l|2007h model predicts V- 
band luminosity- weighted stellar metallicities that are 0.5 
dex below our nw model at all scales, 0.1-0.2 dex higher than 
our vzw model at all scales, and 0.2-0.5 dex below our cw 
model at all masses except the blowout scale. The best agree- 
ment is hence with our fav ored vzw model, i n dicati ng that, 
broadly, the effects of the IKobavashi et al.1 (|2007l ) hyper- 
nova feedback model are similar to our momentum-driven 
outfiows although the hypernova winds may be somewhat 
weaker. 



3.3 The IGM Metallicity 

A fully self-consistent model for galaxy evolution must ac- 
count for the distribution of metals in the IGM as well as 
in galaxies, hence we have also tested our outflow models 
by comparing the predicted and observed IGM metallici- 
ties. Because the impact of our outflow prescriptions on the 
distribution of met als in different ggis phgises has been dis- 
cussed previously (lOppenheimer fc Pavel [20061 : iDave et al.l 
[20071), we will only mention the results here. The mean 
metallicity of the IGM at z = 2 in our favored vzw model 
is roughly [Z/H]=-1.7 l|Dave et al.ll2007l . Figure 1). This is 
roughly twice as high as predict e d by the hypernova feed- 
back model of IKobavashi et al.l l|2007l . Figure 18), imply- 
ing once again that the outfiows in the hypernova feed- 
back model are somewhat weaker. However, by analyzing 
simulated quasar absorption spectra generated alon g sight- 
li nes t hrough the simulation volume, ,Oppenheimer fc Pavel 
(1200^ ) have shown that the predicted abundance of CIV in 
the vzw model is in excellent agreement with observations 
from z — 6 2. Hence our favored vzw model is broadly 
consistent with observations of the distribution of metals 
both within and outside of galaxies at 2 = 2. 



4 ANALYTICAL MODEL FOR THE MZR 

Figure [T] strongly suggests that the observed MZR can be 
used to constrain superwind models, but it gives little phys- 
ical insight into how superwinds impact the trends and evo- 
lutionary behavior of the MZR. In order to probe this ques- 
tion more deeply, we construct a simple one-zone chemi- 
cal enrich ment model s imilar to many others in the litera- 
ture f e . g . ■ I Tinsle vll 1 98"ol ) . The novel aspect is that we will use 
our simulations to calibrate the model inputs, so essentially 
by construction our analytical model will broadly reproduce 
the simulation results, as we show in i]6.3l Owing to its sim- 
plicity, it provides an instructive tool to examine the relative 
importance of various physical effects in driving the MZR, 
which we will do in SjTj 



4.1 Equations of Evolution 

At each timestep, the mass of metals in the ISM Mz in- 
creases owing to inflows and star formation and decreases 
owing to outflows and metals being locked up in long-lived 
stars: 

Mz — ZigmMacc + yA/sFR — ZgMsFn. — 2gA/wind 

= azZgMACC+yMsFB - ZgMsFK{l + Vw) (9) 



8 Finlator & Dave 



Here, 

the ISM, respectively; Macc and Msfr denote the gas ac- 
cretion and star for mation r ates, respectively; y denotes the 
net Oxygen yield l|Tinslevl ll98a l. AiTwind denotes the rate 
at which gas enters the wind; and r;w denotes the ratio 
Afwind/AfsFR- Equation [5] makes the following assumptions: 

• The rate at which new metals are injected into the ISM 
is given by j/Msfr, where the yield j/ is a constant; 

• The outflow rate is proportional to A/sfr; 

• The metallicity of the outflowing gas is equal to the 
mean metallicity of the galaxy's ISM; 

• We assume instantaneous feedback (|Tinslevl ll98Cj ) so 
that the effects of Type la supernovae and delayed mass 
loss are neglected; 

• The mean metallicity of inflowing gas is some fraction 
Qz ^ of the metallicity in the galaxy's ISM; 

• We neglect the difference between the rate of formation 
of all stars and the rate of formation of low-mass stars, i.e. 
we assume that the bulk of stellar mass is in long-lived stars; 
and 

• All metals resulting from star formation and gas accre- 
tion are assumed to be well-mixed within the galaxy's ISM. 

The first four assumptions mimic those made in our simu- 
lations whereas the latter three are made for convenience. 
Also, note that any effects resulting from the enhanced cool- 
ing rates of metal-enriched gas will be taken into account 
implicitly when we tune our star formation efhciencies to 
match the simulations in Section [5.21 

To obtain the metallicity evolution in this model we 
need the evolution of metal mass as a function of gas mass. 
This requires us to know all the parameters in equation [5] 
except the net yield y (we normalize our metallicities by y 
in all of our Figures). We will assume az ~ for the nw 
and cw models and 50% for the vzw model (we will justify 
this assumption in §[7]3}. For 77W in the vzw model, we infer 
a from the time-dependent relation between baryonic mass 
and velocity dispersion that we measure directly from the 
vzw outputs and then plug this into t^w = 300 kms~^/o-. 
For the cw model, we set rjw = (2, 0) for all galaxies with 
masses (below, above) the blowout scale. 

Finally, we need Macc and Msfr as a function of red- 
shift and galaxy mass. In the following sections we describe 
how we calibrate relations for these quantities from our sim- 
ulations. 



4.2 Gas Accretion History 

We infer the time-dependent rate at which galaxies accrete 
fresh gas from their environments directly from the vzw sim- 
ulations. In detail, we trace each galaxy's baryonic growth 
rate backwards in time by searching for its most massive 
progenitor in each simulation snapshot. The total baryonic 
accretion rate is then simply the baryonic growth rate plus 
losses due to outflows: 



643 2 



Macc = A^bar + ??wA/sfr. 



(10) 



In practice this technique yields the accretion rate of gas and 
stars rather than just gas as it is difficult to determine from 
the outputs of our simulations whether new stellar mass 
results from star formation or mergers. However, baryonic 
mass growth is known to be dominated by gas accretion at 
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Figure 3. Mean baryonic accretion rates of the simulated galax- 
ies m the 32/i-lMpc vzw simulation as a function of time (bot- 
tom) and redshift (top). Four different mass bins are shown, with 
the labels referring to the stellar mass at 2 = 0; each history is 
normalized to its maximum accretion rate. Galaxies of all masses 
experience remarkably similar baryonic accretion histories. 



high redshifts (|Keres et al.ll2005l : IGuo fc White! l2007h . and 
the bottom panel of Figure [14] suggests that this is approx- 
imately true in our models as well. 

To estimate outflow losses for the vzw model, we need 
the average r/w for each galaxy as a function of time. We can 
obtain a good estimate by combining the known baryonic 
mass with the time-dependent relation between galaxy ve- 
locity dispersion and baryonic mass, measured directly from 
our simulations, and substituting this back into the relation 
assumed by our vzw model rjw = (300 kma~^ /a). Figure [3] 
plots the resulting accretion histories in 4 bins of stellar 
mass, where we have normalized each one to its maximum 
accretion rate; the mean stellar masses in each bin at z = 2 
are indicated. 

Galaxies of all masses experience remarkably similar ac- 
cretion histories in our simulations, with accretion rates ris- 
ing steadily at early times until a peak around z ~4-3 and 
falling afterwards. The existence of this generic accretion 
history has two interesting implications. First, under the as- 
sumption that Msfr tracks Macc (which we justify below) , 
the rising accretion rates immediately explain why galax- 
ies generically exhibit rising s tar formation histori es at high 
redshifts in our simulations (|Finlator et al.l [20071 ). Second, 
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the existence of a generic accretion history implies that a 
correlation between Msfr and halo mass — and, under rea- 
sonable assumptions, stellar mass— i s exp ected, as found in 
other simulations bv iFinlator et al. (120061) and recently ob- 
served at 2 ~ 0.5 — 2 ( Noeske et al l2007l : E. Daddi, private 
communication) . 

In detail, there are slight differences after z = 2 with 
a suggestion that the most massive galaxies continue to ac- 
crete too much gas at late times. Future work incorporat- 
ing a treatment for AGN feedback is expected to alleviate 
this problem. Additionally, note that the use of a completely 
generic gas accretion history along with the assumption that 
the gas processing rate exactly tracks the gas accretion rate 
does not reproduce the so-called "downsizing" of galaxy 
evolution, where the latter is defined as the tendency for 
more massive galaxies to exhibit lower specific star forma- 
tion rates (,Zhcng ct aj. 2007; Iglcsias-Paraino ct al. 2007), 
unless the ratio of past-averaged to present wind suppression 
factors (1 -I- ?7mlf)/(1 -I- ??w) scales strongly with mass. It is 
more likely star formation is delayed in low-ma ss galaxies by 
an effe ct that our model does not account for dNoeske et al.l 
l2007bl ). 

We averaged over the normalized accretion histories in 
Figure [3] and found that a reasonable approximation to the 
resulting mean accretion history is given by the fitting for- 
mula from ISpringel fc HernquistI l|2003br ): 



Macc{z) cc 



&exp[a(z — Zm)] 
b — a + a exp[b{z — Zm)] 



(11) 



Using a least-squares algorithm we determine the best-fit 
parameters for the vzw model as a = 1.06, b — 1.32, and 
Zm = 3.5. Although this formula was originally designed to 
describe the global star formation history, as iKeres et al.l 
l|2005l ) has pointed out, this closely tracks the gas accretion 
history, so it is not surprising that a good fit is obtained for 
Macc using this form. 

Similar trends are seen in the cw and nw models 
at 2 ^ 2 (unfortunately only our vzw simulations were 
evolved to z = 0). In detail, there are slight differences be- 
tween the gas accretion histories in the three wind mod- 
els. The most prominent of these is that gas accretion 
rates peak at an earlier redshift in the cw model than in 
the vzw or nw models beca use the cw's ene rgetic winds 
heat the IGM too efficiently (|Oppenheimer fc' Dave 2006). 
As a result, the gas cooling time at moderate overdensi- 
ties (~ 10) surpasses a bubble time and the star forma- 
tion ra te density begins to decline at an earlier redshift 
{z ~ 5: ISpring el fc Hernauist"2Q03b') than in the vzw or nw 
models {z ~ 3; Oppcnhcimcr & Dave 2006). We account for 
this early peak by setting Zm = 5 when modeling the impact 
of cw outflows in our analytical model. We will show in Fig- 
ure [12] that using this parameterized fit to the vzw model's 
gas accretion history allows us to approximate the chemical 
evolution of our simulated galaxies in all three wind models, 
hence a more detailed discussion of the impact of outflows 
on gas accretion histories is beyond the scope of the present 
work. The amplitude of a galaxy's gas accretion history is 
determined by a constant multiplicative factor. 
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Figure 4. Least-squares fit to the vzw star formation efficiency 
at z = 2. Because the slope of the trend does not vary with scale, 
we can readily use it to tune our analytical model. 



4.3 Star Formation Rates 

To obtain Msfr, we measure Mspn/Mg in the simulations 
as a function of baryonic mass and redshift, and then use 



MSE 



MsF 



Ma 



-{z,Mb 



(12) 



where the gas mass Mg increases owing to inflows and de- 
creases owing to outflows and star formation according to 



Mg = A/acc - AfsFR ( 1 + w ) • 



(13) 



Through trial and error we determined that it is not possible 
to reproduce the detailed mass-metallicity evolution of our 
simulated galaxies without such a calibration for each wind 
model; indeed, this is a hint as to what governs the MZR. 

Figure 13] shows that the star formation efficiency varies 
with scale in a simple way in the vzw model at 2 = 2; simi- 
lar trends hold for all epochs and wind models. In order to 
tune our analytical model, we fit regression lines of the form 
log(MsFR/Afg) = a -I- 61og(Mbar) to each simulation at a 
range of epochs. Figure [5] gives the resulting fits. In all three 
wind models 6 > 0, indicating that more massive galaxies are 
generically more efficient at converting their gas into stars 
owing to their higher gas densities; this behavior is qualita- 
tively similar to the idea that sta r formation timescal es are 
shorter in more massive galaxies (jNoeske et al.ll2007bl l. The 
offsets a in the wind models are everywhere lower than in the 
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Figure 5. Least-squares fitting parameters to the star formation 
efficiency as a function of baryonic mass in our tliree wind mod- 
els. For eadi model, the Une indicates how fits to the relation 
log(MsFR/Afg) = a + 61og(Mbar) evolve over cosmic time. 



nw model because superwinds decrease galaxies' gas densi- 
ties. The slopes b axe everywhere higher in the wind mod- 
els because wind effectiveness generically scales with mass. 
However, the fact that the cw efficiencies lie very close to the 
nw efficiencies at high redshift and then diverge from them 
suggests that cw winds are relatively ineffective in the low- 
mass galaxies that dominate at early times. The fact that 
the slope and offset of the trend evolve in opposite direc- 
tions regardless of wind model indicates that gas densities 
in low-mass galaxies are more sensitive to changes in envi- 
ronment than in more massive galaxies. It is intriguing that 
the slope of the evolutionary trend changes sign near the 
point z — 0.5, likely a consequence of the universe becoming 
A-dominated. 

Does self-enrichment impact our simulated star forma- 
tion efficiencies? If gas in more massive halos cools more 
efficiently owing to enhanced metallicities then a scale- 
dependent "positive feedback" could obtain between the gas 
cooling and self-enrichment rates; in this case the form of the 
MZR could be inffuenced by the relative effectiveness of this 
positive feedback cycle as a function of mass. Comparing 
Figures [T] and [S] shows that, indeed, the nw model exhibits 
both the highest star formation efficiencies and the highest 
metallicities, in qualitative agreement with this picture. 

However, the low star formation efficiencies in our wind 



models can also be explained by the tendency for gas parti- 
cles to be ejected by winds before their densities (and hence 
star formation rates) grow comparable to the typical den- 
sities in the nw model. Moreover, while the metallicities 
in the nw model scale more steeply than the vzw model, 
the star formation efficiencies scale more steeply in the vzw 
model, in qualitative disagreement with the self-enrichment 
picture. Additionally, we note that galaxy growth at these 
redshift and mass scales i s dominated by cold-mode gas ac- 
cretion (|Keres et al.ll2005l) . with the result that the gas cool- 
ing timescale is much shorter than the dynamical timescale 
irrespective of its metallicity. In Section[7]we will show that, 
in the outflow model that reproduces observations, the MZR 
can be understood entirely in terms of the effects of outflows. 
Hence while a self-enrichment feedback cycle must be oper- 
ating on some level, its effects on the observable MZR are 
Ukely to be weak compared to the effects of outflows. 



5 THE EFFECTS OF WINDS 

In this section we begin our exploration of the impact of 
outflows on the MZR by evaluating how winds affect our 
simulated galaxies; in short, we wish to determine what out- 
flows do to galaxies. In our analytical model we assume that 
the primary parameters through which winds modulate the 
observable MZR are the star formation efficiency, 77W, and 
the gas accretion history. We would like to compare how 
each of these parameters scales with mass in our three wind 
models. Th e gas accretion history has been explored else - 
where (e.g. lOppenheimer fc Davei[200^ : iDave et all l2006al ) , 
and is reasonably well described by equation 1111 Hence we 
first discuss the impact of outflows on r/w and star forma- 
tion efficiency. Next, we show that the fraction of metals 
retained by galaxies does not drive the MZR even though it 
is affected by outflows. Afterwards we compare how outflows 
suppress stellar mass and gas-phase metallicity on a galaxy- 
by-galaxy basis. Finally, we discuss how winds impact the 
trajectories that galaxies follow through the MZR. 



5.1 Mass Loading 

In our simulations, the rate at which material enters the 
wind is given by Msfr?7w. In the cw model, rjw ~ 2 
and Vw ~ 484 kms~^ while in the vzw model r;w ~ 
300 kms~^/(T and Vw ~ 6.7(t. After a wind particle leaves 
the star forming region it interacts with the galaxy's halo 
hydrodynamically. Some gas particles may escape into the 
IGM while others may radiate away their kinetic energy and 
rain back down as "galactic fountains" . Hence it is not obvi- 
ous what fraction of the gas that enters a galaxy's wind will 
actually escape from the galaxy permanently, or whether 
this fraction will preserve the input ryw scaling. Fortunately, 
the connection between star formation, metal enrichment, 
and winds allows us to constrain this quantity. 

If an isolated star-forming galaxy generates a wind with 
a constant mass loading factor rjw then the fraction of metals 
that it retains is given by 



Mz 



= 1 - 



?7w 



(14) 



Mz,formcd 1 + 77W V 

If the masses and metallicities of a galaxy's gas and stellar 
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Figure 6. Mean mass loading factor experienced by galaxies as a 
function of stellar mass. In the nw model only low-mass galaxies 
lose any baryons, while in the cw and vzw models galaxies the 
mass of baryons lost scale as expected. Note how rapidly winds 
become ineffective above the blowout scale in the cw model. 



low masses that is visible in all three of our models. Be- 
cause these effects are small compared to the overall trends 
that relate to the MZR (for example, compare the size of 
the rjiMLF in the nw model to the typical t^mlf in the vzw 
model), they do not affect any of our results. 

In the cw model, below the blowout scale t^mlf ~ 1.5 « 
jTjw, indicating that on average | of wind particles from 
these galaxies leave permanently. Above the blowout scale 

10^°Mq at z ~ 2 - 3) ^MLF declines rapidly because 
these winds thermalize their kinetic energy efficiently and 
return to the source galaxy as galactic fountains. Such a 
phenomenon, if real, would leave clear signatures in the ob- 
served baryonic TuUy-Fisher relation or MZR at this scale. 
As we discuss in ^ the absence of such features in observa- 
tions argues against the cw model. 

In the vzw model, 77mlf oc M^r"'^^, shallower than the 
predicted slope of —1/3 if the stellar mass is a fixed fraction 
of the halo mass. To understand this, in Figure [7] we show 
the stellar fraction /« as a function of Afhaio, and indeed 

1/3 

we see that for the vzw case /, oc Mj^^j^. Taking this into 
account, we find that t^mlf oc M^^^^^ oc 1/a as expected. 
This indicates that outflow processes preserve the assumed 
scalings once the scaling of /» is accounted for as long as the 
galaxy is above the blowout scale. Furthermore, it indicates 
that in the vzw model the fraction of wind particles that es- 
cape the galaxy is roughly constant for all galaxy masses, as 
required by obse rvations of the low -redshift baryonic TuUy- 
Fisher Relation (|Geha et aLlbOOd ). 

It is interesting to note that the r;MLF trend in our vzw 
scenario is qualitatively simi lar to the trend of wind strength 
versus mass in Figure 16 of lKobavashi et al.l l|200'iD . In our 
work, this trend results directly fr om the assumed s caling 
of the instantaneous r/w while in iKobavashT et al.1 (|2007f ) 
it results from their treatment for pressure-driven outflows 
from supernova and hypernova feedback. 



phases are known then Equation [14] can be solved for the 
instantaneous mass-loading factor rjw- In general, rjw varies 
as the galaxy grows so that it is not possible to recover 
rjw from our simulation outputs. However, in this case it 
is still possible to obtain an "effective mass loading factor" 
?7mlf as a function of stellar mass in our simulations from 
Equation [T?] using the known masses and metallicities of the 
simulated galaxies; in this case r^MLF measures the mass loss 
rate averaged over the gala xy's star formatio n history. 

As previously noted bv lDalcantonI l|2006l ). the fact that 
equation [14] is not in general proportional to the effective 
yield indicates that the relative metal contribution to the 
IGM from different galaxies cannot straightforwardly be 
inferred from their effective yields (as attempted by, e.g. 
iBouche e t al. 2007). 

Figure [6] shows t^mlf vs. M, in our various wind models 
at z = 2. This plot is one of the most important ones in this 
paper for understanding the MZR. Looking at the nw model 
first, we find that ~ 10% of the nw galaxies show evidence of 
having lost some of their baryons (t^mlf > 0). These losses 
owe primarily to tidal stripping, although additional scat- 
ter is introduced by discreteness effects and uncertainties in 
the identification of low-mass galaxies within the simulation 
outputs. The result is a slight "fiaring" of the trends toward 



5.2 Star Formation Efficiency 

We now discuss the impact of outfiows on the star formation 
efficiency. We seek to answer two questions: (1) What impact 
do outfiows have on the integrated star formation efficiency 
/, (defined as the fraction of baryons in a halo converted to 
stars as a function of halo mass) at z = 2?; and (2) How do 
outflows suppress /«? The latter question is of interest be- 
cause the expected scalings depend on whether outflows cou- 
ple efficiently with the remaining gas in the galaxy's ISM as 
well as its ambient IGM. If coupling is poor then, for galaxies 
above the blowout scale, the /, scaling should refiect the in- 
trinsic rjw scaling. Alternatively, if coupling is efficient then 
the /* scaling should refiect the condition that the feedback 
energy is comparable to the binding energy of the baryons 
in the halo (e.g., Dekel & Sffk, 1986; Dokcl & Woo 2003). 

Figure]?] (similar to Figure 5 of lDave et al.ll2006al ) shows 
how /* varies with stellar mass in the different wind mod- 
els at 2 = 2. Examining the nw model first, we see that 
in the absence of winds /« climbs steadily with mass below 
A^haio = W^^Mq and then decreases slowly with increasing 
mass above the characteristic minimum mass scale for ha- 
los to be dominated by hot gas rather than cold gas (e.g., 
iBirnboim et al.ll2007l ). This shape ciualitat ively mimics the 
behavior predicted bv ]Dekel fc Wool l]2003l '). However, since 
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Figure 7. Fraction of baryons converted to stars as a function of 
halo mass (bottom axis) and circular velocity (top axis) at 2 = 2. 
The two loci correspond to halos from the 16 and 32/i~^Mpc vol- 
umes that contain more than 128 star particles. The dashed and 
solid lines in the wind models show the range of scalings expected 
from c onsiderations of mas s-loading and energy balance, respec- 
tively llDekel fc Wool |2003| . see text). Below the mmimum mass 
scale for hot-mode gas accretion /, scales with mass even with- 
out outflows. In the presence of outflo ws, ,f« is determined by the 
combi ned effects of mass- loading and lScannapieco fc Broadhursli 
ll200lh suppression. 



there are no winds to couple the feedback energy to the halo 
gas, the qualitative agreement is merely a coincidence. The 
fact that /* is not constant with stellar mass in the nw model 
has the important implication that /, (and hence the MZR) 
is no t governed solely by outflows (Dalcanton 2006). Wc note 
that lTassis et al.l (|2006, Figure 2) and lDe Rossi et al.l (|200^ . 
Figure 2) have observed qualitatively similar behavior in the 
absence of strong outflows. 

What determines the /* scaling in the absence of out- 
flows? If all halos converted their gas into stars with the same 
instantaneous star formation efficiency MsFii/Mg then nei- 
ther Zg nor /* would scale with mass. Evidently Msfr/A^s 
must scale with mass at some point prior to z = 2. In- 
deed, we find that, without outflows, star formation effi- 
ciency and gas density increase strongly with increasing 
mass before the reionization epoch z ^ 6 although both 
trends weaken significantly by 2: = 2. The scaling in star 
formation efficiency at high redshift can be understood in- 
tuitively as a consequence of the fact that more massive 



halos begin collapsing earlier than less massive halos, giv- 
ing them a "head-start" in condensing their gas reservoirs. 
The same effect also manifests itself in a trend for more 
massive galaxies to exhibit older mean stellar ages than less 
massive galaxies in this model, dubbed "natura l downsiz- 
ing" bv lNeistein. van den Bosch, fc Dekel I (|2006h . As a re- 
sult, low-mass galaxies possess lower gas densities, enhanced 
gas fractions and suppressed gas-phase metallicities with re- 
spect to massive galaxies. 

The vzw model qualitatively resembles the nw model 
in /* vs. Mhaio except that it is shifted down by a factor 
of 5-12 with an additional dependence on mass (as shown 
more clearly in Figure |9} . The normalization is lower be- 
cause the vzw model significantly delays star formation (and 
suppresses gas-phase metallicities) at all scales. The fiatten- 
ing behavior in halos above lO^^'^M© obtains because more 
massive halos are generally hot mode-dominated, just as in 
the nw model. 

What determines the slope at masses below this scale? If 
we assume that, as the galaxy grows, the fraction of baryons 
that forms stars at each mass scale is 1/(1 + r/w), then /, is 
simply given by 



/ dM/(l + T7w) 
JdM 



(15) 



This scaling is denoted by the dashed line in the figure, 
where we have made the approximation that the halo ve- 
locity dispersion is given by g = 0.0083(M/Ma )^/^ kms"\ 
Alternatively, if we follow iDekel fc Wol ()2003[ ) and assume 
that stars continue to form until the total energy in outflows 
equals the virial energy of the halo's baryons then we find 
that /. cx M^^'-'il+zY^^ in the vzw model; this scaling is de- 
noted by the solid line in the figure and heis been normalized 
to 1 at the reheating scale at this redshift M = 5.6x10^^Mq. 
At a glance the 77w-driven explanation seems more accurate 
although both theories produce the correct scaling. However, 
noting that the normalizations of these two scaling relations 
are somewhat uncertain owing to the assumptions in our 
spherical collapse estimates, in practice either explanation 
could be valid. Below we show that the /« scaling that is 
expected from energy considerations is not obeyed in the cw 
model, implying that energy in outfiows does not couple effi- 
ciently with infiowing material. For this reason we conclude 
that the accuracy of the /* scaling that is predicted from 
energy considerations is purely a coincidence and that /* is 
dominated by the scaling of the mass loading factor rjw- 

In the cw model /, is relatively fiat between the small- 
est halos containing ^ 128 star particles and the blowout 
scale at a halo mass of lO^^'^M©. Above the blowout scale 
/, climbs slowly. The way that /« varies only slowly with 
mass below the blowout scale and climbs above it echoes 
the cw model's MZR (Figure [!}, indicating a connection 
between the suppressed star formation efficiencies and the 
suppressed metallicities. However, the fact that Figures [T] 
andQlook rather different in the cw model (especially above 
the blowout scale) indicates that the suppressed /» does not 
by itself determine the MZR. Additionally, the fact that the 
MZR and the /* plots look qualitatively different from the 
nw and vzw models indicates that some process specific to 
the cw model is affecting /« and Zg in similar ways. 

What determines how /« scales in the cw model? If /» 
depended only on rjw then we would expect /* = 1/3 at 
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all scales; this is indicated by the dashed line in the Figure. 
The slow dependence of /* on M agrees qualitatively with 
the expected flat trend. However, the normalization is much 
lower than expected, indicating that another process must 
be sup pressing galaxy growt h. Returning to energy consider- 
ations, iDekel fc Wool l|2003l ) used the assumption that stars 
continue to form until the total energy in outflows equals the 
virial energy of the halo's baryons to predict that, in energy- 
driven wind scenarios such as our cw model, /» (and hence 
the metallicity) should increase smoothly with increasing 
mass below the reheating scale of ~ 150 kms~^ and then 
flatten out above it. The range of possible scalings that they 
derive is indicated by the solid lines, where we have normal- 
ized their scalings to /* = 1 at the reheating scale. Their 
predictions were based on the assumption that energy in 
outflows couples efflciently with the baryons in the galaxy's 
ISM and its halo. It is clear from Figure [7] that this key as- 
su mption does no t hold in our cw model. At high masses 
the iDekel fc Wool l|2003l ) model predicts too little suppres- 
sion, indicating that outflows suppress star formation even 
if much of the gas does not escape from the galaxy's halo. 
Furthermore, at low masses their model predicts too much 
suppression, indicating that simulated outflows either escape 
these galaxies without entraining the bulk of the halo gas 
(as would be expected if the outflows are not spherically 
symmetric) and/or their energy is thermalized and radiated 
away. Thus, despite the impressive agreemen t between ob- 
servat ions and the simple model put forth bv lDekel fc Wool 
l|2003t) , simulations suggest that the interaction between out- 
flows and ambient gas is qualitatively different than what 
they assumed. 

We can speculate as to why go o d agre ement with the 
MZR is obtained by iDekel fc Wool (120031 ) . as well as in 
se mi-analytic mod els based on this type of scenario such 
as lDe Lucia et ahl (2004). In these models, they assume that 
winds from galaxies below the reheating scale inject progres- 
sively more energy per baryon into ambient gas, thereby 
unbinding progressively more of it to smaller masses. In 
essence, they force a scaling of 77mlf with Af, that is sim- 
ilar to our vzw model (cf. Figure (Gjl, by assuming that 
j?MLP 3> 77w for small masses owing to efficient energy 
coupling with ambient gas. This physical process is not 
borne out by our three-dimensional numerical simulations, 
in which reheating is mostly irrelevant. It is a minor co- 
incidence that in constant wind models, the blowout and 
reheating scales are fairly similar (equations [T] and [Sjl. 

In summary, a comparison between Figures [T] and [7] in- 
dicates that (1) in the presence of outflows, the scaling of 
/* (or, equivalently, gas fraction) depends more heavily on 
the mass-loading factor ryw than on energetic considerations; 
and (2) /, , while a key driver, does not by itself determine 
the MZR. Some other factor must be important in determin- 
ing the basic shape of the MZR. We will argue below that 
it is primarily gas accretion, although gas stripping due to 
winds from neighboring galaxies can play a role in some sit- 
uations, as we will discuss in In ij I5.3l and 15.41 

5.3 Metal Retention 

The most popular interpretation of the observed MZR is 
that low-mass galaxies exhibit low metallicity because they 
drive a larger fraction of their metals into the IGM. In order 
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Figure 8. Fraction of metal mass retained fz,ret as a function 
of M, . The two loci correspond to galaxies from the 16 and 
32h~'^Mpc volumes. Without winds galaxies tend to retain their 
metals, whereas in the presence of winds galaxies can lose up to 
50% of their metals. These trends are far too weak to account for 
the observed MZR. 

to determine whether this idea explains the observable MZR 
in our simulations, we plot in FigureOthe fraction of metals 
retained fz,iat as a function of stellar mass in each of our 
wind models. If metal loss dominates the form of the ob- 
servable MZR then we expect the slope and scatter of /z,ret 
as a function of M, to mimic the MZR. For simplicity we 
discuss only the 16/i~^Mpc boxes at z = 2 here while noting 
that similar results hold for other scales and epochs. 

In the nw model galaxies retain all of their metals on 
average (as expected in the absence of outflows) although 
~ 0.06 dex of scatter is introduced by dynamical disruption 
and uncertainties in group idenflcation. This contrasts with 
the strong scaling {Zg cx M^'^*) and larger scatter (0.11 
dex) in the nw MZR (Figure [!}. Moreover, a few low-mass 
galaxies have lost up to 80% of their metals through interac- 
tions, yet we flnd that even these galaxies exhibit no depar- 
ture from the mean gas-phase metallicities for their stellar 
masses. Hence while they have ejected significant baryons, 
those baryons were enriched at the same level as the baryons 
that remained in the galaxies. 

Turning to the wind models, at a fiducial stellar mass 
of 10^° M0 the cw and vzw models retain roughly 40% and 
70% of their metals sd z — 2, respectively. However, the 
mean gas-phase metallicity of the cw galaxies is roughly 70% 
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higher than it is in the vzw galaxies. Comparing the scalings 
reveals similar inconsistencies: In the vzw model /z.rot oc 
M° while the MZR scales as Zg oc M°'^\ In the cw model, 
/z.ret oc M-°-°'^ with 0.09 dex of scatter, while its MZR 
scales as Zg oc M°"^ (note that this is in the opposite sense 
as the /z,ret scaling) with 0.19 dex of scatter. 

The poor correspondence between the /z.rct scalings 
and the MZRs for all three models indicates that galaxy 
metallicities are not primarily driven by metal loss. The in- 
consistency between the relative gas-phase metallicities and 
/z.rot values of the vzw and cw models suggests that the de- 
tailed way in which metals are distributed in different bary- 
onic phases must be taken into account. Note that we do not 
claim that fz, ret cannot trace the MZR m principle; Figure[8] 
only shows that it does not do so in general. Indeed, there 
is no rigorous reason why it should. Observationally, it is 
possible to test this if the metallicities and masses in the 
stellar and gas phases as well as the net metal yield y are 
known: if Zg oc /z.rot then it should also be proportional 
to Zt{y — Mg / M,)~^ at all scales. However, guided by our 
own simulations, we will argue in 33 that the MZR is dom- 
inated by the scaling of the star formation efficiency in the 
absence of winds, and by a competition between the rates 
of enrichment and dilution in the presence of winds. 



5.4 Suppression of K'h and Zg 

Another way to highlight the differences between the wind 
models is to compare their effects on individual galaxies. Be- 
cause all of our simulations were run with the same initial 
conditions, we can readily do this by matching the positions 
of individual galaxies between the simulations. Figure [9] dis- 
plays the ratios of stellar mass and metallicity in the wind 
models versus the no wind model at 2 = 2. Note that we 
have directly verified that excluding from Figure |9] those 
galaxies whose nw metallicities suggest t^mlf > does not 
impact the slope, normalization, or scatter of the inferred 
suppression factors for either wind model. 

In the cw model one might at first expect the ra- 
tios M«,cw/Mt,nw and .Zg,cw/^9,nw to be « 1/(1 -I- ?7w) for 
galaxies below the blowout scale as long as outflows cou- 
ple inefficiently with the ambient gas (i.e. if rjw = 2 then 
2/3 of the baryons should be ejected in winds). Figure [§] 
showed that cw outflows are highly effective for all galax- 
ies below the blowout scale, in that their effective mass 
loading factor is, on average, equal to 3/4 of the true as- 
sumed r;w = 2, hence we might expect Af*,cw/Mt,nw to 
be roughly 1/(1 + 2.5) = 0.4. Turning to Figure]^ we see 
that Mt_cw/M«,nw lies below this value even at the blowout 
scale, consistent with the existence of an extra source of 
suppression implied by Figure [7] Above the blowout scale 
the cw ratios climb with tight scatter, reaching unity at 
M, K, lO^-^M©. 

Below the blowout scale M«,cw/M,,nw varies slowly with 
decreasing mass and shows considerable scatter. At the low- 



est masses (Af, < IQ Mq) Af,,cw/Mt,n„ and Zg,cvj / Zg^-nvj 
actually climb with decreasing mass despite the fact that 
jjMLF is constant in this range. We have verified that these 
trends continue to lower masses at higher mass resolution, 
hence the behavior is not an artifact of numerical resolution. 
Moreover, it clearly conflicts with the nai've picture above. In 
order to understand it, we must consider how winds might 
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Figure 9. Ratio of stellar mass and gas-phase metallicity in the 
wind models versus the no- wind model at 2 = 2. The two loci 
correspond to galaxies from the 16 and 32/i~^Mpc volumes. In 
the cw model the effect of wind suppression does not vary strongly 
below the blowout scale whereas it decreases rapidly above it; in 
the vzw model wind suppression is less effective at higher masses 
and the scaling is quite smooth. 



affect galaxy growth in the full context of structure forma- 
tion rather than as isolated systems. 

First, there is the possibility that an early genera- 
tion of galaxies pre-enriches the IGM, giving rise to a 
"metallicity floor" in the observable MZR. Such a mini- 
mum has been inferred from quasar absorption line sys- 
tems (e.g., 'Songailal l200ll) as well as in ou r simulations (e.g. 
iDave et al. 2007). In IPave etHI ()2006bl ') we hypothesized 
that widespread pre-enrichm ent of the IGM in the cw model 
(jOppenheimer fc Dave|[2006l ) may be responsible for its low- 
mass behavior. One way to quantify pre-enrichmcnt is to 
measure how many metals were created in one galaxy, ex- 
pelled, and reaccreted onto another. We measured the frac- 
tion of galaxies' metals at 2; = 2 that originated in another 
galaxy and flnd that it is roughly 5% and 15% in the cw and 
vzw models, respectively. While this is not insignificant, it 
is not sufficient to account for the scatter in the cw MZR. 
Furthermore, the fact that the vzw model has more pre- 
enrichment but shows much less scatter suggests that this 
process does not contribute significantly to the MZR scat- 
ter. Hence the pre-enrichment hypothesis is unlikely to be 
correct. 
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Second, in a full hierarchical context there are two ef- 
fects that reduce the tendency of winds to drain galaxies' 
gas reservoirs, which may lead to increased metallicities. 
Whereas in a simple closed-box scenario all gas is equally 
eligible to form stars, in a fully self-consistent model galax- 
ies' gas reservoirs possess density and pressure gradients and 
star formation is concentrated in relatively small regions. 
Gas that is ejected from these areas can in principle be re- 
placed on a dynamical time by infall from less dense areas 
owing to the loss of pressure support. A second effect is 
that gas entering a wind can potentially thermalize its en- 
ergy before escaping the galaxy's halo and fall back down 
onto the galaxy in a galactic fountain, effectively increasing 
the galaxy's gas accretion rate. Figure |6] does not indicate 
the presence of scale-dependent galactic fountains below the 
blowout scale. On the other hand, the fact that M«,cw/M«,nw 
increases at low masses suggests that gas that escapes the 
galaxy in a wind is indeed being rapidly replaced, boosting 
Mt,cw/Af*,nw above what would naively be expected. 

Finally, galact ic winds can affect neighb o ring galax- 
ies. In particular, IScannapieco fc BroadhurstI (120011 ) pro- 
posed that halos with masses 1O®~^°M0 suffer stripping 
of their baryons owing to winds from lower-mass neigh- 
bors. In their model, once galaxies form they drive spherical 
winds whose energies are comparable to the wind energies 
in our cw model. The winds from galaxies in low-mass ha- 
los (< l(f Mq) then strip baryons from intermediate-mass 
halos (10^~^°Mq) that are not yet virialized without af- 
fecting more tightly-bound massive halos (> 10^° Mq), with 
the result that galaxy formation in intermediate-mass halos 
is suppressed (note that this effect is not related to tidal 
stripping). Although our winds are not spherical, when av- 
eraging over a large sample of galaxies the effect should still 
be noticeabl e. It is expected to be weak at 2 ^ 6 (e.g.. 
Figure 7 of iDave et al.l l2006al ) and to grow most notice- 
able dur ing the heyday of galaxy forma, tion z ^ 3 (Fig- 
ure 8 of IScannapieco fc BroadhurstI I2OO1I ') . Indeed, the ra- 
tios M*,cw/M,,nw and .Zg,cw/^9,nw show significantly less 
suppression at 2 = 6 (not shown) than at z = 2. More 
interestingly, they also show less scatter at higher red- 
shifts, as expected for an effect whose strength depends 
on environment rather than on galaxies' intrinsic proper- 
ties. In summary, it is likely that stellar mass growth in the 
cw model divides into three mass regimes: At low masses 
[Mt/MQ < 10*) galaxy growth is suppressed only by out- 
flows and hence M*,cw/M»,nw — > l/(l+ryw). At intermediate 
masses (10* < Mt/M© < 10^'') galaxy growth is dominated 
by a collaboration between IScannapieco fc BroadhursH 
l|200ll ) suppression and losses to winds, where the rel- 
ative contribution of each effect varies nontrivially with 
scale and e poch. Because it only works in a particular 
mass regime, IScannapieco fc BroadhurstI (I2OOII ) suppression 
hence gives rise to a local minimum in Mt,cw/M,,nw. Fi- 
nally, at high masses (M^/M© > 10^°) galaxy growth is 
again dominated by outflows, where the effects of outflows 
weaken rapidly with increasing mass owing to the rapidly 
declining »7mlf (Figure [S]). 

In the vzw model wind suppression of metallicity and 
stellar mass both scale smoothly, with relatively little scat- 
ter, and without evidence for a preferred scale. We have 
already shown (Section 15. 2p that stellar masses in the vzw 
model can be attributed largely to the scaling of 77W. Given 



that fz.Tct does not scale strongly with mass in this model, 
it is also not surprising that Zg^vzw/^g.nw scales in the same 
way. 

Does baryonic stripping occur in our vzw simulation? 
We can predict whether it should be stronger or weaker in 
the vzw versus the cw model by estimating the ratio of their 
momentum generation rates: 

Pvzw ,vzw?7W,vzw V'w,vzw (16) 

yew ,cw 

From Figure [T] /,,vzw ~ 0.1 and /,,cw ~ 0.2. From be- 
fore, we know ^7w,vzw — 

300 kms-^/a, Vw, VZW 

6.7a, 

2, and Vw.cw = 484 kms ^, hence we flnd that 
^vzw/fcw ~ 1 — the effect should be roughly as strong in 
the vzw model as in the cw model. By subtracting the 
mean trend from the vzw MZR and inspecting the resid- 
ual, we have found that vzw galaxies less massive than 
M* — 10^'^ Mq do, in fact, exhibit signatures of baryon strip- 
ping, seen as a change of slope and a slightly increased scat- 
ter in the MZR below this scale. However, the effect is much 
less noticeable than in the cw model because it is small com- 
pared to the effects of the strong outflows. Applying this 
idea to interpret observations, the tight sca tter and smooth 
scaling in the observ ed low-redshift MZR (iLee et al ] |2006l : 
iTremonti et al.|[2004r ) indicate that either galactic winds do 
not carry enough momentum to suppress growth in neigh- 
boring halos appreciably or the effects of outflows must be 
strong enough to erase any evidence of baryon stripping; 
both possibilities clearly conflict with the behavior of our 
cw model, but are satisfied in our vzw case. 

Comparing the top and bottom panels of Figure [5] 
yields insight into the extent to which /, governs the MZR. 
The metallicity suppressions show considerably more scat- 
ter than the mass suppressions. Additionally, the gas-phase 
metallicities are systematically about 50% less suppressed 
than stellar masses in the cw models, independent of M«. 
(They are roughly equally suppressed in the vzw model.) 
Both of these observations are inconsistent with the idea 
that /* solely governs the MZR; evidently the effects of 
winds on gas accretion rates and star formation efficiencies 
also play a role. 

In summary. Figure [5] indicates that outflows tend to 
suppress both Af» and Zg. The trends and the levels of 
scatter indicate that the amounts of suppression cannot, 
in general, be predicted from simple scaling relations. We 
have highlighted the fact that, in both flgures, the scatter 
in the cw model seems to be large below the blowout scale 
and small above it, whereas the scatter i n the vzw model 
is comparable to --^ 0.1 dex at all masses. iLee et all (|2006l ) 
recently reported that the la scatter in the observed rela- 
tion is roughly 0.1 dex from lO^-lO^'^M© at low redshift 
and argued that this implied a less energetic form of metal- 
enhanced mass loss than blowouts. Our results tend to sup- 
port this view. However, we also find that, in order to avoid 
the la rge scatter introduced by IScannapieco fc BroadhurstI 
(l200lt ) suppression, winds must either transport signifi- 
cantly less momentum out of halos than occurs in the sim- 
plest energy-driven wind models, or else they must invoke 
large mass loading factors 77W > 2 at intermediate masses. 
Momentum-driven winds naturally satisfies that require- 
ment while our constant wind model does not. 
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Figure 10. Mean enrichment histories of the simulated galaxies 
in the three wind models starting from 2: = 6 — > 2. Blue and red 
solid lines denote evolutionary trends from the 16- and 32/i~^Mpc 
boxes, respectively, while the green dashed line gives the slope of 
a closed-box with baryonic massSxlO^M© {Zg/y = — ln/i(M«)). 
The magenta curves indicate the mean observable trends at 2 = 2. 
Galaxies evolve roughly parallel to the observable trend in the vzw 
model and more steeply in the cw and nw models. 



6 EVOLUTION OF THE MZR 

In the previous sections we explored how the effects of winds 
scale with mass by studying galaxies at a single epoch. In 
order to understand how outflows affect the evolution of the 
MZR, we now discuss how the metallicity evolves through 
cosmic time in our three models, as a function of galaxy mass 
and gas fraction. To do so, we trace the evolution of each 
simulated galaxy at z = 2 by searching for its most massive 
progenitor in each simulation output, thereby compiling its 
mass and metallicity history. Because enrichment histories of 
individual galaxies are highly stochastic, we bin the histories 
by stellar mass in order to show the typical evolution as a 
function of stellar mass. 



6.1 Evolution of M, ~ Zg 

Figure [To] compares the mean enrichment histories from the 
different models from 2 = 6 — * 2. As expected, galaxies gen- 
erally increase in both mass and metallicity as they evolve. 
However, the slope of the evolution is in general neither con- 
stant nor parallel to the observable trend owing to the fact 



that gas accretion rates, star formation efficiencies, and wind 
properties vary with scale and time. In particular, galaxies 
do not generally evolve as closed boxes (green dashed line) 
although the nw model comes quite close even at high red- 
shift. Instead, their evolution is more shallow owing to the 
fact that outflows expel gas that is enriched compared to 
inflows from the IGM. At lower redshifts where strong out- 
flows are rare and accretion rates are low, the evolution is 
expected to more closely resem ble the closed - box s cenario. 
This can be seen in Figure 1 of iBrooks e\~^ (|2006| ). where 
their simulated galaxies evolve from z = 2^0 with a slope 
dlog(Zg)/dlog(M.) « 1. Additionallv. lSavaglio et al.l (|2005l ) 
have shown that closed-box scenarios can account for the 
observed evolution from z — 0.7 0.1. 

There is a slight upturn in the evolutionary slope from 
2 = 3-^2 that can be seen in aU three models. This fea- 
ture occurs because gas accretion rates are declining dur- 
ing this interval owing to the increasing gas cooling times 
at intermediate overdensities; in other words, because gas 
accretion grows decreasingly effective at diluting gas reser- 
voirs, star formation grows increasingly effective at enriching 
them. The fact that the rate at which the normalization of 
the MZR changes during this interval varies among our wind 
models illustrates how the time evolution of the MZR, in ad- 
dition to its normalization, shape, and scatter, is a testable 
prediction of the wind model. 

Turning to the individual models, galaxies in the nw 
model evolve nearly as steeply as a closed box in all but 
the highest mass bin. The relatively steep evolution in the 
low-mass bins reflects the fact that these galaxies' metallic- 
ities are not significantly diluted by inflows owing to their 
large gas fractions. By contrast, the shallower evolution in 
the highest mass bin indicates that, owing to their high gas 
densities and star formation efficiencies, these galaxies pos- 
sess somewhat lower gas fractions with the result that in- 
fiows dilute their metallicities more readily. We will show in 
Figure[TT]that the same effects can be seen in these galaxies' 
effective yields. 

In the vzw model galaxies enrich their gas reservoirs 
somewhat more slowly. This evolution can intuitively be un- 
derstood as a consequence of the tendency for vzw galaxies 
to expel a large fraction of their metals. However, we can 
also understand it via our analytical model as follows: In 
the vzw model, the effective mass loading factor scales with 
stellar mass as t^mlf oc M^"'^^ (Figure[S|. If galaxies' metal- 
licities remain near equilibrium (Section 17. ip and if the gas 
processing rate is approximately equal to the gas accretion 
rate Msfr(1+?/w) ~ ^acc, then the vzw MZR should scale 
with mass as Zg ~ y/{l + rp/^). For rjw S> 1 (typical of small 
galaxies at high-2) we therefore expect Zg cx M^'^^ — exactly 
the scaling seen in Figure [10] 

The evolution of the MZR in the cw model is more 
complex than in the other two models owing to the nontrivial 
environment-dependent interactions between the outflowing 
winds and the inflowing gas. The relevant points are that the 
basic shape of the constant wind MZR does not evolve with 
redshift while its normalization increases somewhat more 
rapidly than that of the vzw but not as rapidly as the nw 
model. 
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Figure 11. Mean enrichment histories of the simulated galaxies 
in the three wind models starting from z = 6 ^ 2. Dashed green 
lines show where galaxies fall that have ycff/y = 0.5, 1.0, 2.0. Blue 
and red solid lines denote evolutionary trends from the 16- and 
32h~^Mpc boxes, respectively. The nw galaxies involve much like 
closed boxes while the vzw and cw galaxies show suppressed ef- 
fective yields. 



6.2 Evolution of Effective Yields 

A more direct way to investigate the extent to which our 
simulated galaxies depart from closed-box evolution is to 
plot the metallicity versus gas fraction, Figure [11] The axes 
in this figure are chosen so that galaxies with a constant 
effective yield j/cff = Zg/ ~ ln(/i) (where p denotes the mass 
fraction of baryons in the gas phase) evolve along straight 
lines. The dashed green lines show the evolution for (from 
bottom to top) Vcs/y = 0.5, 1.0, 2.0. Individual galaxies gen- 
erally evolve from lower left to upper right in this space. 

The nw galaxies remain quite close to the closed-box 
curve ycs/y = 1-0 as suggested by the steep evolution in 
Figure [TO] The effective yield ne ver exceeds the c losed-box 
case, as required by Theorem 3 of lEdmundsl l|l990l ). It drops 
below the closed-box value most strongly at the redshift cor- 
responding to the peak accretion rates z ~ Z because the 
dilution of the metallicities at this epoch is overcompensat- 
ing for the increase in gas fraction (this happens irrespective 
of the wind model and is therefore a robustly predicted — if 
difflcult to confirm — consequence of the global gas accretion 
peak at z ~ 3). After the accretion rates begin subsiding. 



2/off/y returns to the closed-box value relatively q uickly ow- 
ing to continued star formation l|DalcantorJl2006l l . 

By contrast, galaxies in the vzw model tend to evolve 
from higher to lower j/cff/y during this period owing to the 
combined effects of dilution and outflows. In this model ac- 
cretion rates fall off somewhat more slowly following the 
z = 3 peak than in the other mod els, pushing the peak 
of st ar formation to lower redshift (jOppenheimer fc Pavel 
l2006f ) and delaying the recovery of the effective yields. The 
delayed falloff in accretion owes to galactic fountain effects 
that occur because outflowing gas does not escape the halo 
~ 20 — 50% of the time, on average. The gas is then retained 
in a puffy distribution, owing to small wind speeds that do 
not drive gas far from the galaxy, and becomes available for 
reaccretion on a time scale smaller than a Hubble time. 

Galaxies in the cw model show a gradual evolution from 
nearly closed-box yields Vcs/y ~ 0.8 at z = 6 to lower values 
Vas/y ~ 0.3 at z = 3 before rebounding to ycs/y ~ 0.5 at 
2 = 2. The relatively weak dependence of ycs/y on M» at 
all redshifts recalls the flat MZR in Figure [T] and owes to 
the fact that most of the galaxies in this figure lie below 
the blowout scale and are hence roughly equally affected by 
winds. 

The bottom-right plot in Figure [TT] compares how ^eff 
varies with M* at 2 = 2 in our three models vers us observa- 
tions of UV-selected galaxies at 2 2 (|Erb et al.i 2006). The 
nw model shows no trend at low masses and declines with in- 
creasing mass for masses above 10^°Mq. The decline occurs 
because gas fractions decline with increasing mass: Galax- 
ies more massive than 10^^ Mq have gas fractions /i 20% 
in the nw mod el so that unenri ched inflows are able reduce 
T/off efficiently (|Dalcantonll2006l ). The yetr behavior is qual- 
itatively consistent with the observed trend although the 
normalization is ~ 50% too high. 

In the vzw model, j/cft increases with increasing mass 
at low masses and fiattens out around 10^^ Mq. This behav- 
ior at face value conflicts with observations at 2 ~ 2. On 
the other hand, it bears a striking resembl ance to the ob- 
served trend at low redsh ift (Garnctt 2002: FlYemonti et al.l 
l2004l : |Pflvugin et al.ll2004 ). Ahhough we do not show it here, 
the overall shape of the j/eff trend in the vzw model does 
not vary with redshift. An aggressive interpretation of Fig- 
ure [TT] would be that winds at high redshifts must differ 
qualitatively from winds at low redshifts, with the former 
relatively more effective in massive galaxies than the latter; 
in other words, at high redshift the effect of winds on j/eff 
should preserve the nw scaling while at low redshi ft the vzw 
model is more realistic. However, considering that lErb et al.l 
(2006) were forced to infer gas masses indirectly rather than 
measuring them, and given that their measurements span a 
much smaller dynamic range than low-redshift observations, 
we prefer not to draw any firm conclusions from Figure [TT] 
Future measurements that trace high-redshift gas masses 
more directly (e.g. with ALMA) will constrain the effective 
yields of high-redshift galaxies more directly, and will pro- 
vide a key test of the momentum-driven wind scalings. 

The vzw trend at low masses is expected despite their 
higher gas fractions (not shown) because t^w increases to 
lower masses (Figure [SJ . In terms of our analytical model, 
we now anticipate the conclusions of Section [7] by assuming 
that Msfr(1 -I- Tjw) = AfACC and dividing Equation 1201 by 
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yln(l//i) to find the equilibrium condition: 

J/cff,cq _ 1 



y 



(17) 



(l + r,w)ln(i) 

As mass increases, r]w decreases (Figure [Gjl, hence yoff.cq 



mcreases. 



The flattening behavior in i/cff,cq around A/, = 
1O^'^~^°M0 is more interesting. Recall that in our vzw 
model the fraction of baryons lost to winds does not vary 
strongly with mass (Section 15. 3|l . Evidently, a flattening in 
1/cff does not necessarily indicate a scale at which superwind 
feedback becomes effective at removing a galaxy's metals, 
an interpretation we discuss further below. Nor does it in- 
dicate that j/eff has reached the true yield yes/y ~ 1- In- 
stead, it seems to be a coincidence resulting from a com- 
petition between dilution owing to accretion and the out- 
flows' suppression of star formation. Indeed, in more mas- 
sive galaxies (M* > W^^ Mq) ycs begins to increase with 
mass again owing presumably to the fact that gas accretion 
rat es are falling in this ra nge. While the highest mass bin in 
the lTremonti et al.l (j2004l ) data tantalizingly suggests agree- 
ment with this behavior, confirmation will have to await the 
arrival of larger samples of massive star-forming galaxies at 
2 ~ 1 (since at 2 = massive galaxies are generally too 
gas- poor for t heir g as-phase abundances to be measured). 

iGarnetd (|2002l ) also found a fiattening of y^a below 
~ 100 kms~^ at low redshift, but his interpretation was 
that low-mass galaxies retain a smaller fraction of their 
metals owing to the onset of winds in less mas sive galax- 
ies, in conflict with our interpretation. However, IPalcantonl 
recently used the same data to show that the ob- 
served trend in j/efp can be obtained without assuming that 
winds are more effective in low-mass galaxies by removing 
the approximation that the gas fraction is constant. In this 
view, i/off is suppressed in low-mass galaxies because their 
low star-formation efficiencies prevent their i/eff from recov- 
ering quickly from inflow episodes. Our findings tend to sup- 
port the latter view. 

In the cw model, ydt does not vary with mass below the 
blowout scale and increases with increasing mass above it, 
emphasizing the idea that the effectiveness of winds below 
the blowout scale does not seem to follow the trends that 
would be expected from simple scaling arguments. 

Effective yields can also give insight into another impor- 
tant question in the study of high-redshift galaxies, namely 
whether the buildup of stellar mass occurs in a predom- 
inantly smooth or episodic fashion. At redshifts z < 1, 
the tight observed correlation between stellar mass and 
AfspR argues in favor of a predominantly smooth mech- 
anism for relatively massive galaxies 10^° < M,/Mq < 
10^^ ijNoeske et al.1 l2007h . At z ~ 2, clustering measure- 
ments indicate that Lyman-break and submillimcter galax- 
ies po ssess duty cycles of ^ 1 ([Adclbcrgcr et al. 2005a) and 
~ 0.1 (|Bouche et al.ll2005l ). respectively, suggesting that Ly- 
man Break galaxies form stars smoothly while submillimeter 
galaxies are more bursty. Similarly, recent clustering mea- 
surements indicate that Lyman-a emitters at z ~ 4.5 pos- 
sess duty cycles of ~ 10 %, again hinting at relatively bursty 
star formation histories l|Kovac et al.ll2007l ). 

Effective yields provide another way to test this behav- 
ior, because in the absence of inflows and outflows, galaxies' 
effective yields quickly recover to the true yield. Hence if the 



timescale for the yas/y to recover to the true yield is short 
compared to the duty c ycle then some m echanism must be 
actively suppressing it I DalcantonI l200d ) . Quantitatively, if 
we define the ratio of the effective yield to the true yield 
Xy = ycff/y then, in the absence of inflows and outflows, 
Xy varies with time according to 



dXy 



MsFH. 



dt A41n(i: 



■{l-Xy 



(18) 



This equation s hows that the equilibrium solution Xy = 1 
is a stable one (|Koppen fc Edmunds! [l999h and that depar- 
tures from equilibrium disappear with an e-folding timescale 
given by Mg ln(l/^)/AfsFR. Applying this timescale to UV- 
selected galaxies at 2 ~ 2, where by computing weighted 
means over bins 2-6 in lErb et al.l ^006 we determine 
(< AfsFR >,< Mg >,< /X >) « (29 A/Qyr~\2.1 x 
10^" Mo, 0.38), we find that the timescale for yes/y to re- 
turn to the closed-box yield is 700 Myr. This is consider- 
ably shorter than the gas consumption time of 1.2 Gyr. This 
short timescale together with the suppressed observed effec- 
tive yields j/cff < 0.01 imply that star formation in these 
galaxies cannot be episodic in nature, consistent with the 
large inferred duty cycles. In our simulations, this timescale 
is less than 50% of the Hubble time at all epochs in both 
wind models, indicating that our predicted effective yields 
reflect star formation and gas accretion processes that are 
predominantly smooth rather than episodic in nature. 

In summary, at high redshift {z = 6 —> 2) galaxies 
tend to evolve from high to low j/eff as accretion rates in- 
crease and are hence expected to "rebound" at lower red- 
shifts z < 2 where gas accretion rates decrease. Superwind 
feedback generically suppresses y^B in our galaxies as ex- 
pected. The resulting Af* — ycB trend depends on the ways 
in which winds affect how inflows, outflows, and gas frac- 
tions scale with mass. The cw trend conflicts with high- and 
low-redshift observations while the vzw trend agrees quali- 
tatively with low- but not all high-redshift observations. The 
fact that the vzw model produces better agreement with the 
directly-measured high-2 MZR and \ow-z Af» — ycB trends 
suggests that our vzw model is closer to reality. In both cases 
star formation occurs via predominantly smooth rather than 
bursty processes. It will be interesting to see whether the 
vzw model's conflict with the observed trend in ycB at high- 
z is alleviated once direct measurements of gas densities at 
high redshifts become available. 



6.3 Analytic Model MZR Evolution 

In Section [7] we will synthesize the insights gained from our 
analytical model to understand the origin of the MZR in our 
simulations. In order to justify this analogy, in this section 
we show that our analytical model broadly reproduces the 
evolution of the simulated MZR. 

Figure [12] compares how galaxies evolve through the 
MZR in our simulations (thin lines) versus our analytical 
model (thick lines) for three different mass bins. The analyt- 
ical model succeeds in recovering most of the qualitative fea- 
tures of the models as well as the differences between them. 
For example, in each model the observable MZR has nonzero 
slope. This is only achieved if the star formation efficiency 
or the riw scales with mass (Section 17. 2|) . The agreement 
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Figure 12. A comparison of the mean evolution of galaxies 
through the MZR in the wl6 (thin blue) and w32 (thin red) mod- 
els versus our analytical model. The analytical model reproduces 
the qualitative differences between the wind models although it 
does not reproduce the behavior of the individual models in detail. 
In particular, the analytical model yields a nonzero slope at low 
masses in the cw model, in clear conflict with the hydrodynamic 
simulations. 



the gas-phase metallicity. The upturn at late times is not 
conspicuous in any of our analytical models although it can 
readily be reproduced by forcing gas accretion rates to drop 
more precipitously at 2 < 3 than we have done. This is es- 
pecially cleax in the cw model, which is generally the most 
difficult model to reproduce owing to the fact that its behav- 
ior does not scale smoothly with galaxy mass. Our analytical 
cw model inevitably yields a nonzero slope in the observable 
MZR, in conflict with the simulations. This can be alleviated 
by allowing the star formation efliciencies MsFn/Mg to vary 
more slowly at low masses than at high masses (Figure [13)) . 
However, further fine-tuning of the star-formation efficien- 
cies and accretion histories in our analytical model would 
yield little insight. The important point to take away from 
Figure [12] is that our analytical model captures the essential 
ingredients that determine how galaxies evolve through the 
MZR in our fully three-dimensional simulations. In the next 
Sections, we will therefore apply our analytical model to de- 
termine the conditions that drive the form of the observable 
MZR. 



7 UNDERSTANDING THE 

MASS-METALLICITY RELATION 

In this section we use the intuition gained in the past several 
sections along with the simulation results in order to piece 
together a comprehensive understanding of what drives the 
MZR's form and evolution. We have already demonstrated 
that our vzw simulation produces good agreement with the 
slope, normalization, and scatter of the observed MZR (Fig- 
ure [TJ and that our analytical model provides an accept- 
able description of how galaxies evolve in our full simulation 
(Figure [12]). Hence we begin by showing how our analytical 
model can account for the amplitude, slope, and scatter of 
the simulated — and, by implication, the observed — MZR. 



between the simulated and analytical trends results because 
we directly tuned our analytical star formation efficiencies 
to match Figure [5] We have verified that our analytic model 
reproduces the simulated galaxies' effective yields as well 
(Figure [TT|) although we do not discuss it here. 

The qualitative differences in evolution between our dif- 
ferent outflow models are well-reproduced. For example, nw 
galaxies evolve the most from jz = 6 — > 2 while cw galaxies 
evolve the least because the nw galaxies retain their entire 
gas reservoirs while the accretion histories in the cw model 
peak at an earlier time. The vzw galaxies evolve the most 
shallowly owing to their high r/w's while the nw galaxies 
evolve the most steeply. The result of these differences is 
that the vzw galaxies are the least enriched at 2: = 2 while 
the nw galaxies are the most enriched. 

While our analytical model reproduces many of the 
gross features of the simulations as well as the differences 
between the different wind models, it does a poorer job of 
reproducing the individual wind models in detail. For ex- 
ample, in the nw and vzw models metallicities at z = 2 are 
« 0.1 dex too low. Such offsets can easily result from mi- 
nor inconsistencies in accretion histories or the deflnition of 



7.1 Implications of the Model 

7,1.1 Normalization and Scaling 

Combining our original analytical model, equation [9] and 
the evolution of the gas mass, equation 1131 it is straightfor- 
ward to show that galaxies evolve through the MZR with a 
slope given by 



AZg ^ J_ / Macc 
AM, Mg VMspr 



which is equal to zero if 



y 



MsFR 



Macc(1 - ctz) 



(19) 



(20) 



This possible balance between the influences of star for- 
mation and infall ha s been identified previous ly (e.g., 
Tinslev fc Larsonlll978l : iKoppen fc Edmundj|l999D . Winds 
enter into the determination of Zg^eq indirectly by modulat- 
ing the rate at which a galaxy depletes its gas reservoir as 
well as the relative enrichment of the satellite galaxies that 
it accretes. If a galaxy processes its gas into stars and winds 
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at the gas accretion rat^B, then Msfr(1 + ?7w) = Macc, 
which in turn yields the equilibrium gas-phase metallicity 
•^s.cq = t/(l + r;w)~^(l — az)~^ ■ In what follows we will 
show that in our wind models this is approximately true. 



7.1.2 Scatter 

The ratio of a galaxy's metallicity to its equilibrium metal- 
licity Xz = Zg/Zg^cq evolves with time according to 



dt 



ACC 



Ma 



(1 - Xz) - yXz 



Ms 



Ms 



Macc \ 
Macc / 



(21) 



The second term in Equation[2T]is small except during short- 
lived interactions, hence we may neglect it. In this case, 
Equation [21] i mplies that the equilibrium solution Xz = 1 
is a stable one (|Koppen fc EdmundsllT999l ) and that depar- 
tures from equilibrium disappear on a timescale given by 
Mg/ Macc, or the timescale for the gas reservoir to be di- 
luted by a factor of 2. If this timescale is shorter than the 
timescale over which perturbations to a galaxy's metallic- 
ity occur then galaxies' gas-phase metallicities recover from 
perturbations quickly, suppressing scatter in the observable 
MZR. 



7.2 Normalization and Scaling Without Outflows 

We begin our discussion of the origin of the MZR with the 
no-wind scenario in order to develop some intuition about 
how hierarchical structure formation impacts the MZR. In 
Figure [1] we showed that its MZR has approximately the 
correct slope, but an amplitude that is ~ 0.5 dex too high. 
We further showed in Figure |S] that no- wind galaxies retain 
most of their gas, as expected without outflows. Hence the 
only physical effect left that can cause a slope in the MZR 
is the star formation efficiency. 

If galaxies converted their gas into stars at precisely 
the gas accretion rate then the slope of the MZR in the 
nw model would be z ero. Th i s can be seen from the "infall 
model" formalism of iLarsonl (| 19721 ) . which tells us that if 
the gas mass is constant and Msfr = Macc then the gas 
metallicities evolve as: 



u = fl - 1, 



(22) 



(Note that this can be obtained from Equation [19] by sub- 
stituting Msfr/Macc = 1.) In this model, for sufficiently 
small gas masses the MZR would be flat because the gas frac- 
tion /i would rapidly shrink to zero for all galaxies. Equiva- 
lently. Equation [21] tells us that we would expect all galaxies 
to approach Zg^cq — y as long as Mg /Macc were signifi- 
cantly less than the bubble time. The MZR clearly is not 
flat in the nw model. Therefore, guided by our discussion of 
Figure [7] we now ask whether the scaling of the star forma- 
tion efficiency MsFn/Mg can account for the scaling of the 
no-wind MZR. 

We have calculated the z = 2 MZR using our analyti- 
cal model with three different prescriptions for the SFR. In 
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Figure 13. Predicted MZR at 2 = 2 in the nw model versus 
several analytical models. The solid black curve gives the mass- 
weighted gas-phase metallicity in the nw simulation. The red 
dashed and blue dotted lines are computed by assuming that gas 
collapses into stars in 1 and 10 dynamical times, respectively (see 
text); the long-dashed magenta line is computed by tuning star 
formation efficiencies to match the nw model. Green dot-dashed 
curves denote infall models (Equation I22I I assuming constant gas 
masses of 10*, 10^, and 10^^ Mq from left to right. The slope of 
the MZR in the absence of winds is dominated by the scaling of 
the star formation efficiency. 



the first, we compute the star formation rate by assuming 
that gas cond enses into s t ars in 10 dynamical times; this is 
similar to the iKennicuttI (|l998l ) relation. The gas densities 
are obtained from the baryonic masses by combining the 
Virial T heorem with the l ow-redshift baryonic Tully-Fisher 
relation (jCeha et al.ll2006l ) and assuming that galaxies are 5 
times as dense as their host halos (the value chosen for this 
ratio affects the amplitude but not the trend of our results) . 
In the second, we assume that gas condenses into stars in 
one dynamical time. In the third, we tune the star formation 
efficiencies Mspn/Mg to reproduce the efficiencies in the nw 
simulations. 

Figure [13] compares the resulting MZRs at z = 2 with 
the trend from the fully three-dimensional nw model using 
mass- weighted gas-phase metallicities (solid black curve). 
Additionally, we have plotted the evolution of Zg versus M» 
for three representative infall models (Equation 1221 green 
dot-dashed curves) corresponding to constant gas masses of 
10**, lO'^, and 10^°Mq from left to right. The trend from 
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the 10-dynamical time prescription has a negative slope, 
indicating that more massive galaxies possess larger gas 
masses (and larger gas fractions) than less massive galax- 
ies and hence cannot enrich their gas reservoirs as efltec- 
tively; in terms of Equation 1191 since Msfk/ Macc declines 
with increasing M*, Zg does as well. The trend from the 
one dynamical time model is flat, as expected for a scenario 
with such efficient star formation; indeed, for this model 
MsFR ~ Macc so that Equation [22] accurately predicts 
Zg/y — 1 for all masses since the gas fractions are negligi- 
ble. By contrast, the trend from the model in which we have 
tuned the star formation efficiencies exhibits the desired pos- 
itive slope. This trend is in qualitative agreement with the 
numerical trend although the most massive nw galaxies ex- 
hibit metallicities above the yield. The high metallicities at 
the massive end are likely a consequence of accreting pre- 
enriched gas; in terms of Equation 1191 az > 0. 

The simulated trend is shallower than expected from 
Equation [22] because the assumption Msfr ~ Macc is 
violated weakly in the absence of outflows; in particular, 
Msfr/Macc drops from 1.0 at W^Mq to 0.9 at lO^M© 
in the analytical (magenta) curve while drops from 0.6 
to 0.2 over the same interval. In short, the reason the nw 
model has the correct MZR slope is that the star formation 
efficiency in a hierarchical structure formation scenario nat- 
urally yields the desired differential with galaxy mass (cf. 
§[S3}. Nevertheless, the excessive amplitude (cf. Figure [TJ 
suggests that metals must be preferentially removed from 
these galaxies. He nce outflows are n ecessary to obtain the 
correct amplitude ijKobavashi et all [20071. but are not re- 
quired to obtain the su ppressio n of metallicity i n low-mass 
galaxies ijPe Rossi et al.ii2006i ; iTassis et al.ll2006l ). 

7.3 Normalization and Scaling With Outflows 

Next we consider the impact of winds. From Figure [5] we 
see that outflows suppress both the stellar mass and the 
metallicity. In order to lower the MZR amplitude, outflows 
must lower Zg more than M^'^, which is the observed slope 
of the MZR. 

Looking first at M* = IO^^A/q, which is around L, at 
2 = 2, we see from Figure [T] that both the cw and vzw 
models produce roughly the correct MZR amplitude (recall 
that the observed metallicities and the yield each introduce 
~ 0.3 dex of uncertainty). Under the assumption that the 
MZR is governed primarily by the mass loading factor r/w 
via the equilibrium conditions (cf. 17.1. ip . this is expected 
because galaxies in both models experience effective mass 
loading factors in the range ?;mlf = 1 — 1-5 at this scale 
(cf. Figure [61 . Intuitively, ?7mlf determines the level of sup- 
pression of gas enrichment and stellar mass growth because 
higher mass-loading factors lead to lower gas densities (cf. 
Figure [5} and to more gas being ejected in winds rather 
rather than converted into stars. This leads to a tendency 
for galaxies with similar ?7mlf to possess similar properties 
regardless of wind model. 

Above this mass scale, constant wind outflows cannot 
escape halos, causing t^mlf to drop rapidly. Hence cw metal- 
licites grow much more rapidly with M* than observed. Be- 
low this scale, metallicities reach a minimum in the range 
M* = W^'^ Mq and then "rebound" slightly to lower masses. 
If this trend continues to z = (as it does from z « 6 — » 2), 



then this would be in gross conflict with observations. The 
cause of this scaling likely owes to the flat trend of jjmlf with 
mass at low masses, with the slight rebound at low masses 
due to baryon stripping f rom nearby galaxies' outflows as 
discussed before (see also IScannapieco fc Broadhurst|[200ll , 
Figure 12). 

In the vzw case, galaxies are always below the blowout 
scale (which increases linearly with mass). In this regime, 
the mass loading factor governs the MZR. In order for the 
scaling of Zg oc M°'^ to hold, 77W must be roughly propor- 
tional to M^"'^, which is satisfled in the vzw case. There 
is also slight evidence for baryon stripping, but it is highly 
subdominant compared to the high r^MLF in the mass range 
where stripping is effective. 

In summary, winds suppress galaxy masses and metal- 
licities primarily by modulating the relative rates at which 
gas reservoirs are enriched and diluted and secondarily by 
stripping baryons from neighboring galaxies. In the vzw 
model the latter effect is small compared to the former while 
in the cw model both are significant. We emphasize that the 
normalization and scaling of the MZR are not determined 
by the total fraction of metals that galaxies retain because 
the rates of gas accretion and star formation are too rapid 
for the gas reservoirs to retain any memory of this quan- 
tity. If baryon stripping is negligible compared to the effects 
of outflows then the most important parameter is the ef- 
fective mass loading factor fjuLF', if the wind speed exceeds 
the escape velocity then fjuLF tx -qw and the scaling of r/w 
dominates the MZR. 

7.4 Normalization and Scaling: The Equilibrium 
Metallicity 

We now illustrate our main result that the mass loading fac- 
tor governs the MZR below the blowout scale by applying 
the analytical relations derived in § 17.11 to the vzw simu- 
lation. In particular, we show that our simulated galaxies' 
metallicities do indeed track the equilibrium metallicity. 

For each simulated galaxy in the vzw model, we have 
used our progenitor lists to track how the ratio of its 
gas phase metallicity to its equilibrium metallicity {Xz = 
Zg / {yMsFB./ Macc)', cf. Section r7.1.2p varies with mass and 
time during the interval 6 < 2 < 2. If galaxies' gas-phase 
metallicities closely track equilibrium and if infalling mate- 
rial is unenriched then we expect Xz « 1 at all masses and 
redshifts. 

The result is shown in the top panel of Figure [14] Com- 
paring the least and most massive galaxy bins, we find 
that whereas the actual gas metallicities of W^^'^Mq and 
10* '' M0 galaxies differ by ^ 0.6 dex at 2 = 2 (Figure [l|, 
their Xz's only differ by f» 0.2 dex. Moreover, the spread 
in Xz is even tighter before z = 2. This would not be ex- 
pected in the absence of an equilibrium condition. The fact 
that the spread in these ratios is tighter than the spread 
in galaxies' actual metallicities indicates that at all times 
galaxies' metallicities are tightly constrained by a balance 
between enrichment from star formation and dilution from 
inflows. The fact that the ratios are offset from zero implies 
that the mean metallicity of inflowing gas is more than 10% 
of the mean metallicity in the galaxy's ISM and reflects the 
widespread presence of galactic fountains in the vzw model; 
the ratio az can be read from the y-axis on the right side 
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Figure 14. (Top) Mean ratio of actual to equilibrium gas-pliase 
metallicity in four bins of M* in the vzw model. The different 
curves correspond to galaxies in different mass bins; the legend 
gives the stellar masses at z = 2. The vertical dashed line denotes 
2 = 2. (Bottom) Mean ratio of the gas processing to gas accretion 
rates. The mean metallicity of infalling material grows with both 
galaxy mass and time and is roughly az = 50% of the mean ISM 
metallicity at z = 2 (right y-axis on top panel). Gas processing 
rates generally lie within 20% of the equilibrium values at all 
times. 

of the top panel. The increase in az with cosmic time re- 
flects the growing relative contribution of galactic fountain 
gas with respect to pristine ISM gas: by z = 2, az ~ 50%, 
indicating that roughly 50% of the infalling gas is galactic 
fountain material. The increase in az with mass reflects the 
fact that the rate at which pristine IGM gas accretes onto 
the galaxies declines with increasi ng halo mass owing to their 
increasing hot gas fractions (e.g.. iBirnboim et al.ll2007l ). 

The top panel of Figure [TJ] also supports our view that 
pre-enrichment of gas that falls onto galaxies is not signif- 
icant in the vzw simulation, which can be understood as 
follows: If, at a given redshift, gas that is being accreted 
were homogeneously pre-enriched to a certain level, this pre- 
enrichment would provide a relatively larger boost to the 
metallicities of low-mass galaxies than to massive galaxies. 
As a result, the gas-phase metallicities of low-mass galax- 
ies would lie farther above their expected equilibrium value 
given unenriched infall, and the normalization of the Xz 
trend for low-mass galaxies would be boosted systematically 
above the Xz trend for more massive galaxies. In fact, all 



but the most massive galaxies display roughly the same ra- 
tio Xz down to 2 = 2, hence a homogeneous pre-enrichment 
is not significant in this model. 

In Sections [5] and [7] we proposed that the rate at which 
gas is processed into stars and winds tracks the rate at which 
it is accreted, Msfr(1 -|- »7w) — A/acc- This idea is central 
to the current work as it allows us to demonstrate that the 
slope and normalization of the MZR depend almost entirely 
on the scaling of r/w- Additionally, a systematic imbalance 
between the rates of geis accretion and gas processing could 
in principle mimic the effects of nonzero az, leading to an 
incorrect interpretation of the top panel of Figure 1141 For 
this reason, we show the ratio of the gas processing to the gas 
accretion rates in our vzw model in the bottom panel of Fig- 
ure [I4l This figure indicates that galaxies generally process 
their gas at the same rate as they accrete it, justifying our 
assumptior0that Msfr/Macc = 1/(1+^) and supporting 
our view that the increase in Xz with time in the top panel 
results from galactic fountains rather than a mismatch be- 
tween gas accretion and gas processing rates. Interestingly, 
this plot also indicates that, in the vzw model, mergers do 
not contribute significantly to the buildup of galaxies' stel- 
lar populations even in our most massive bin. Because we 
measure baryonic accretion rates rather than strict gas ac- 
cretion rates (the latter being difficult to infer from our sim- 
ulations), a tendency for galaxies to accrete a significant 
fraction of their baryons as ready-formed stars would show 
up as a tendency for the accretion rate to exceed the gas 
processing rate. There is some evidence that this does occur 
in the most massive galaxies in the vzw model, i.e. that dry 
mergers are more prevalent at high masses. However, on av- 
erage galaxies at z ^ 2 do not accrete more than ~ 10% of 
their baryons in th e form of already- formed stars (see also 
iGuo fc Whitll200'iD . 

In summary, the top panel of Figure [TJ] verifies that the 
slope and normalization of the observable MZR are domi- 
nated by the equilibrium condition in Equation l20l while the 
bottom panel verifies that galaxies process newly accreted 
gas into stars and winds at roughly the gas accretion rate. 
In other words, at all masses and redshifts, metallicities are 
dominated by an equilibrium between the rates of enrich- 
ment and dilution while the enrichment rate is dominated 
by an equilibrium between the rates of gas accretion and gas 
processing. These equilibrium rates are governed primarily 
by the mass loading factor, hence the scaling of the mass 
loading factor directly determines the scaling of the MZR. 



7.5 Scatter 

In our analytical model, scatter in the MZR occurs because 
a perturbation, such as an accretion or merger event, dis- 
places a galaxy's metallicity from its equilibrium value. The 
timescale to return to equilibrium is given by the gas dilution 
time Mg / Macc- Perturbations to a galaxy's metallicity are 
expected to occur on timescales no shorter than the dynam- 
ical time, tvir — 2.5Gyr/(l-|-z)^/^. Hence if the dilution time 
is shorter than the dynamical time, then we expect pertur- 



^ Note that this is a ge neralizat i on of the Ansatz MgpR = M\qq 
that was introduced bv lLarsonI lll972l) . 
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Figure 15. The ratio of the dilution time Mg/ M\cc to the 
dynamical time as a function of mass and time in our various 
models. The galaxy masses at 2 = 2 are indicated in the legend 
and the vertical dashed lines indicate 2 = 2. Comparison with 
Figure[T]indicates that mass scales with short dilution times tend 
to show small scatter in the MZR whereas mass scales with long 
dilution times tend to show large scatter. 



dex) in the corresponding MZR: Because perturbations to 
these galaxies' metallicities disappear relatively slowly, they 
spend a relatively large amount of time out of equilibrium. 
Galaxies in the vzw model possess significantly shorter dilu- 
tion times than the nw galaxies or the low-mass cw galaxies, 
consistent with the tight scatter (« 0.08 dex) in the vzw 
MZR at all scales. In the cw model low-mass galaxies pos- 
sess long dilution times and large scatter whereas massive 
galaxies possess short dilution times and small scatter, as ex- 
pected. The most massive cw galaxies possess the shortest 
dilution times of any of our simulated galaxies at z — 2, con- 
sistent with the extremely tight scatter (« 0.06 dex) in the 
cw MZR above the blowout scale. The dilution times for low- 
mass cw galaxies are shorter than they are for the low-mass 
nw galaxies even though the scatter in the MZR is largest for 
the low-mass cw galaxies. The extra scatter at low masses 
owes to these galaxies having had their baryons stripped 
by winds from neighboring galaxies. Because the amount 
of stripping varies with environment, it effectively creates a 
range of equilibrium metallicities for each halo mass. 

In summary, galaxies tend to process gas into stars 
and winds at roughly the gas accretion rate. This tendency 
leads to the existence of an equilibrium gas-phase metallicity 
Zg,eq = J/Msfr/Macc, which encodes information regard- 
ing both a galaxy's present conditions (via the accretion rate 
and wind properties) and its star formation history (via the 
current gas content, which determines A/sfr). Corrections 
to this zeroth-order equilibrium result from galaxies accret- 
ing pre-enriched gas or ready-formed stars as well as effects 
that depend on environment; the first two of these effects 
should not increase the scatter in the MZR while the last 
one should. Metallicities are expected to lie close to Zg,cq 
as long as dilution times are short compared to a dynamical 
time; in our momentum-driven wind model that achieves the 
best agreement with the observed z = 2 MZR, this holds for 
all scales and epochs. Our interpretation of the origin of the 
MZR scatter also explains the no-wind and constant wind 
cases. Together with its relative simplicity, we believe this 
makes our interpretation reasonably compelling. 



bations to disappear rapidly and for the scatter in the MZR 
to be tight. 

Figure fTSl shows how the ratio of the dilution time to the 
dynamical time varies with scale and time in our models. In 
the nw models the dilution times scale inversely with mass at 
early times from roughly l.Stvir in the most massive galaxies 
to roughly IStvir in the lowest-mass galaxies at z ^^2-3. In 
the cw model the dilution time scales more slowly with mass 
although the trend for low-mass galaxies to dilute their gas 
reservoirs more slowly than massive galaxies is preserved 
despite the cw winds. Generally, galaxies less massive than 
10^" M0 have dilution times longer than tvir at z = 2 while 
more massive galaxies have dilution times that are shorter 
than tvii- Finally, in the vzw model the dilution time does 
not exceed fvir at any mass scale. More interestingly, in this 
model it is the low-mass galaxies whose dilution times are 
the shortest, a direct consequence of our assumption that 
the mass loading factor scales with the inverse of the velocity 
dispersion. 

Turning to the question of the scatter in the MZR (Fig- 
ure [l]), the relatively long dilution times in the nw model at 
z — 2 are consistent with the relatively large scatter (~ 0.11 



8 SUMMARY 

In this paper we have compared the observed MZR of 
star-forming galaxies at z ~ 2 with predictions from 
cosmological hydrodynamic simulations that incorporate 
three different models for galactic outfiows: No out- 
fiows, a "constant w i nd" m odel that emulates the energy- 
driven iDekel fc Silk! i 19861 ) scenario, and a "momentum- 
driven wind" mode l that reproduces z > 2 I GM metal- 
licity observations (lOppenheimer fc Davell2006l ). We have 
shown that the momentum-driven wind model produces the 
best agreement with the slope, normalization, and scatter 
of the observed MZR. We have constructed a simple an- 
alytical model that qualitatively reproduces the behavior 
of our simulations, and used it to identify the processes 
that drive galaxies' metallicities. Our work shows that the 
slope, normalization, and scatter of the MZR as well as its 
evolution with time all constitute constraints on the be- 
havior of outflows. In particular, our simulations strongly 
disfavor any constant wind scenario, and explain why our 
momentum-driven wind model produces reasonable agree- 
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ment with available constraints. Our main conclusions are 
summarized as follows: 

• Outflows are required in order to bring the simulated 
and observed MZRs into agreement at 2 ~ 2. Without out- 
flows, enrichment proceeds too rapidly relative to dilution 
with the result that galaxy metallicities are 2-3 x higher 
than observed. 

• The MZR is governed by an evolving equilibrium be- 
tween the enrichment rate owing to star formation and 
the dilution rate owing to gas accretion. This results in 
an equilibnum metallicity for any given galaxy, given by 
Zg,eci = TjMsfk/ Macc, where y is the yield, Msfr is the 
star formation rate, and Macc is the gas accretion rate. 

• Outflows affect ^g,cq, and hence the MZR, by limit- 
ing the gas supply for star formation. For a given mass 
loading factor ryw, and assuming (as our simulations pre- 
dict) infalling gas has a negligible metallicity, then ^g,cq = 
y/{l + Vw)- 

• Wind speeds (Vw) affect the MZR by governing how 
much outflowing mass actually escapes the halo. This re- 
sults in an effective mass loading parameter j7MLF (Fig- 
ure |6]), which is similar to (but slightly less than) r/w so 
long as winds are fast enough to escape a galaxy's halo, 
and drops rapidly towards zero for galaxies whose winds are 
slower than the escape speed. Our simulations' metallicities 
are hence well described by Zg^cq ~ J//(l + Vmlf), where 
?7mlf ~ rjw below the blowout scale, and t^mlf ~ above 
it. 

• The reheating scale, which is the mass below which the 
outflow energy input is sufBcient to unbind all the gas, does 
not play a significant role in determining t^mlf. This is evi- 
dent because in our cw run the fraction of baryons converted 
to stars (/«) does not vary with halo mass in the way ex- 
pected under the assumption of efficient energetic coupling 
of outflows with ambient gas (Figure [7| . Physically, this is 
because in our simulations outflows tend to blow holes in 
surrounding gas rather than heat it. 

• The observed slope and amplitude of the MZR there- 
fore constrain how t^mlf and Vw varies with M, . Our 
momentum-driven wind model obtains the observed relation 
Zg{M:,) oc M°-^ by having rjw ocl/ a oc M^j^^ oc M."^/^ 
and by having outflow speeds always above the escape ve- 
locity (so r^MLF rpfif). The latter constraint requires some 
positive mass dependence of outflow speeds on galaxy mass, 
which we assumed to be Vw oc a, but any dependence where 
galaxy masses are always below the blowout scale would suf- 
fice. 

• Our no-wind scenario also produces a MZR with 
roughly the correct slope, though in detail it is too shallow. 
In the absence of outflows, this owes to a mass dependence 
of Mt on Mhaio such that low-mass galaxies have a lower 
fraction of baryons in stars. Hence in principle it is possible 
to match the observed MZR without having t^mlf vary with 
Af, . Ifowever, the required scaling of M« with Mhaio does 
not occur naturally in our simulations with outflows. 

• Our constant wind scenario fails to even qualitatively 
match the observed MZR. The existence of a blowout scale 
at ~ IO^^Mq produces a marked feature in the MZR, below 
which ?7mlf ~ ?)w ^constant, and above which ?7mlf goes 
rapidly to zero and hence the MZR rises quickly towards 
Zg.cq = V- Such a feature is generically expected across the 



blowout scale. The absence of such a feature in the observed 
z « MZR from M* ^ 10^ - 10^^ M© argues against a 
blowout scale in that mass range, thereby ruling out any 
reasonable constant wind speed scenario. 

• The scatter in the MZR is governed primarily by the 
dilution time td = Mg / Macc compared to the dynamical 
time tvir. If dilution times are short compared to a dynam- 
ical time, then perturbations from Zg^^q have time to equi- 
librate, thereby suppressing scatter. The small scatter seen 
in the MZR argues for td/ivir ^ 1 across the full range of 
observed masses. Our momentum-driven wind model satis- 
fies this non-trivial constraint, whereas our other models do 
not (Figure (TS)). 

• Another physical effect that plays a secondary role 
in governing the MZR is that outflows carry significant 
amounts of momentum that can strip baryons from neigh- 
boring halos. This increases scatter by causing equilibrium 
metallicities to depend on environment as well as mass. 
Hence the tight scatter in the observed MZR suggests that 
either galactic winds do not carry significant amounts of 
momentum out of galaxies, or outfiows must be sufficiently 
mass-loaded to "drown out" the effects of baryonic strip- 
ping. In our constant wind case below the blowout scale, 
neither are true, and the scatter increases significantly. In 
our momentum-driven wind case, the latter is generally true. 

• Outflows and inflows cause galaxies to evolve more shal- 
lowly than closed-box models at early times, with the vzw 
galaxies evolving most shallowly of all (Figure [TTJ. Effective 
yields are expected to be ~ 0.01 at z = 2 for both our wind 
models. However, the detailed scaling of the vzw model's 
2/ofr suggests better agreement with the well-constrained low- 
redshift observations. It is worth noting that the effective 
yield is only reflective of the recent history of gas and metal 
accretion over a dilution time, hence it cannot be used to 
infer long-term accretion histories. 

According to our analysis, an outflow model that will 
successfully reproduce the observed MZR must satisfy three 
main conditions: 

(i) 77W oc(slope of MZR)~^ when ryw ^ 1; 

(ii) Vw must scale with mass such that all galaxies are 
below the blowout scale (so that t^mlf ~ ?7w); 

(iii) Dilution times must be short compared to dynamical 
times in order to maintain a small MZR scatter at all masses. 

These criteria show that the MZR mostly constrains the 
mass loading factor, with weaker constraints on outflow 
speeds and gas accretion rates. It is interesting that our 
momentum-driven wind scenario naturally satisfles these re- 
quirements (along with secondary requirements such as the 
subdominance of baryon stripping). Although other wind 
models could conceivably be postulated that also satisfy 
these requirements, it is compelling that this same model 
also satisfles IGM metallicity constraints, and broadly agrees 
with available direct measurements of outflow parameters at 
high and low redshift. In any case, other wind models will 
likely need to satisfy the above criteria in order to match 
the observed MZR, demonstrating that the MZR provides 
strong constraints on outflow properties. 

The rather dramatic failure of our constant wind sce- 
nario is suprising in light of the app arent success of the 
simple analytical models presented by iDekel fc Silkl (|l986f) 
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and lDekel fc Wool (|2003l ). The root difference traces back to 
those works assuming that feedback suppresses star forma- 
tion by efficiently coupling supernova energy with baryons in 
halos, while our three-dimensional simulations produce in- 
efficient coupling with a propensity for winds to blow holes 
in surrounding gas. This is partly a result of the way we im- 
plement winds in our simulations by turning off hydro forces 
for some distance; however, in practice that distance is much 
smaller than the halo size and hence interactions with halo 
gas can (and do) still occur. Regardless, the existence of a 
strong feature at the blowout scale seems an unavoidable 
consequence in a constant wind scenario, and is in direct 
conflict with the observed unbroken MZR power law over 
four orders of magnitude in stellar mass. Furthermore, the 
increased scatter below the blowout scale predicted by such 
a scenario is not seen. Hence we strongly disfavor this ex- 
planation of the MZR. 

The slow turnover in the 2 ~ MZR at M* > 1O^°-^M0 
cannot be studied directly in our simulations owing to a lack 
of sufficient dynamic range, along with the fact that our sim- 
ulations were not evolved to z = Q. However, it arises natu- 
rally in our scenario when the mass loading factor becomes 
^ 1, which yields ^g.cq ^ y ^constant. In principle, it could 
also arise if galaxies with M, > 10^"'^ M© are above the 
blowout scale (which would also make ?7ml f — > 0); indeed, 
this i s the conventional interpretation (e.g. iTremonti et al.l 
|2004| ). However, this would imply the existence of a blowout 
scale at that mass, which as we have argued above causes 
other features in the MZR that contradict observations. 
Hence we suggest that this mass scale does not reflect a 
characteristic wind speed, but rather a characteristic scale 
of the mass loading factor, namely the galaxy mass where 
the mass loading factor is roughly unity. 

Our findings agree bro adly w i th tho se from the higher- 
resolution study of Brooks et al.l l)2006l l. though there are 
some differences in interpretation. In their work, they de- 
termined that winds affect the MZR of low-mass galaxies in 
the following sense: When they compared gas-phase metal- 
licities at z = with the mean metallicity of all gas that 
had ever belonged to the galaxies, they found no systematic 
offset. This is expected if the mean metallicity in an out- 
flow equals the mean metallicity in the galaxy's ISM. On 
the other hand, by comparing simulations with and with- 
out winds they found that winds suppress gas densities and 
hence star formation efficiencies, which in turn shapes the 
observable MZR. Hence they deduced that star formation 
efficiency is a key driver of the MZR. We also find that 
more massive galaxies have more efficient star formation in 
Figure [S] and Figure [13] shows that this is important for 
establishing the MZR, at least in the no-wind case. 

However, in our wind models we find that the star for- 
mation efficiency doesn't by itself determine the MZR, be- 
cause the trends in Figure [7] don't mimic those of the MZR. 
Instead, the MZR's trend is the mass scaling of Zg^eq, which 
is set by ho w the accretion rate cort ipares with the star for- 
mation rate l|Tinslev fc LarsonI 19781 : eqn. l20|) : in our models, 
this is similar across all masses and close to unity at z; ~ 2, 
when the dependence on r/w is taken into account, as shown 
in the bottom panel of Figure 1141 Hence the mass depen- 
dence in .Zg,eq arises mainly from the mass dependence in 
j?MLF. In this way, the trends in r?MLF (Figure[6| are directly 
reflected in the MZR. 



This interpretation c an be compared with the results 
of iKobavashi et all \2m% . who observed a tight correlation 
between stellar metallicity and the mass fraction of met- 
als retained by the galaxies in their models at all redshifts 
(their Figure 16d) and concluded that higher stellar metal- 
licities result directly from a lower mass fraction of metals 
ejected. Our models obey a similar correlation, hence our 
conclusions should be consistent with theirs. In fact, the 
tendency of gas-phase metallicities to track an equilibrium 
value combined with a tendency for 77MLF to decline and 
Zg,eti to grow as galaxies grow requires that the fraction 
of metals retained and the mean stellar m e tallici ty must 
grow together, as found by IKobavashi et al.l l|2007l ). More- 
over, given that the gas-phase MZR shows little scatter at all 
redshifts, the stellar MZR should also show little scatter (as 
also noted by IKobavashi et al.l 120071 ) even though it is not 
directly governed by an equilibrium condition analogous to 
Eauation l2UI In this way, Equation[5D]likely governs the gas- 
phase MZR directly and the stellar MZR indirectly in both 
sets of models. Indeed, it is likely that an analogous equi- 
librium condition governs the MZR of any galaxy evolution 
model that incorporates a treatment for strong (ryiviLF ^ 1) 
outflows. 

Despite considerable progress over the last decade, the 
use of metallicities and effective yields to constrain galaxy 
evolution constitutes a field that is in its infancy from both 
theoretical and observational perspectives. Our simulations' 
implementation of outflows, while being state-of-the-art for 
cosmological simulations, is still crude. For instance, we cur- 
rently assume enrichment only from Type II supernovae, we 
do not shut off winds in galaxies with low star formation 
rate surface densities ((Heckman , 2003i ) , and we use the local 
potential as a proxy for galaxy mass. All of these simplifi- 
cations are probably not fatal at 2 = 2, but by 2 = they 
likely are; this (in addition to computer time constraints) is 
the main reason we have not attempted to extend our simu- 
lation analysis to z = 0. We are working towards incorporat- 
ing metals from Type la supernovae and stellar mass loss, 
improved wind criteria, and direct galaxy identification into 
our code (Oppenheimer et al., in preparation). Our prelim- 
inary results indicate that none of these issues significantly 
impact the predicted z = 2 MZR. 

Another aspect for future exploration is different scal- 
ings of the wind model. For instance, our constant wind sce- 
nario is only one possible implementation of energy-driven 
outflows. More sophisticated versions that allow 77W and the 
wind speed to vary could improve the agreement between 
the observed and simulated luminosity functions at the faint 
end while yielding agreement with the observed MZR. It is 
by no means clear that momentum-driven winds, as we have 
implemented them, are the only vi able alternative. Indeed , 
it is for this reason intriguing that IKobavashi et al.] (|2007l ) 
have obtained reasonable agreement with the observed MZR 
using a treatment for pressitre-driven outflows that result in 
a qualitatively similar scaling of mass-loading factor versus 
mass as well as the predicted MZR. The intuition gained 
from these results can hopefully guide us (and others) to- 
wards understanding how alternative outflow models may 
fare prior to running expensive simulations. 

On the observational side, galactic outflows are still rel- 
atively poorly constrained despite impressive advances over 
the past decade. More detailed measurements of the mass 
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loading factors and wind speeds across a large dynamic 
range would be helpful in indicating whether momentum- 
driven or energy-driven winds are likely to dominate. More 
importantly, despite heroic observational efforts there re- 
main relatively few constraints on galaxies' mctallicitics and 
gas fractions at high redshift. Upcoming metallicity mea- 
surements made with multi-object infrared spectrographs 
such as FLAMINGOS-2 as well as direct gas mass measure- 
ments made with IRAM and ALMA will prove crucial in fi- 
nally allowing us to apply these metrics to the high-redshift 
Universe. 

Despite its simplicity, the fact that our model explains 
the detailed shapes of both of our wind models' MZRs leads 
us to believe that it captures most of the essential physics. 
Our scenario invokes two parameters, the equilibrium metal- 
licity and the dilution time, neither of which can be directly 
measured because they depend on the gas accretion rate. 
Instead, they must bo constrained indirectly through obser- 
vations of how galaxy properties such as SFR, gas mass, 
and metallicity vary with mass and epoch. We look forward 
to undertaking such comparisons to observations in the fu- 
ture, and are hopeful that they will shed further light on the 
critical problem of understanding galactic outflows. 
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